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ABSTRACT 

Aims. Our goal is to probe the populations of obscured and unobscured AGN investigating their optical-IR and X- 
ray properties as a function of X-ray flux, luminosity and redshift within a hard X-ray selected sample with wide 
multiwavelength coverage. 

Methods. We selected a sample of 136 X-ray sources detected at a significance of > 3er in the 2 — 10 keV band 
(-F2-10 ^ 10" 14 erg cm" 2 s" 1 ) in a ~ 1 deg area in the XMM Medium Deep Survey (XMDS). The XMDS area is 
covered with optical photometry from the VVDS and CFHTLS surveys and infrared Spitzer data from the SWIRE 
survey. Based on the X-ray luminosity and X-ray to optical ratio, 132 sources are likely AGN, of which 122 have 
unambiguous optical - IR identification. The observed optical and IR spectral energy distributions of all identified 
sources are fitted with AGN/galaxy templates in order to classify them and compute photometric redshifts. X-ray 
spectral analysis is performed individually for sources with a sufficient number of counts and using a stacking technique 
for subsamples of sources at different flux levels. Hardness ratios are used to estimate X-ray absorption in individual 
weak sources. 

Results. 70% of the AGN are fitted by a type 2 AGN or a star forming galaxy template. We group them together in a 
single class of "optically obscured" AGN. These have "red" optical colors and in about 60% of cases show significant 
X-ray absorption (Nh > 10 22 cm -2 ). Sources with SEDs typical of type 1 AGN have "blue" optical colors and exhibit 
X-ray absorption in about 30% of cases. The stacked X-ray spectrum of obscured AGN is flatter than that of type 1 
AGN and has an average spectral slope of T = 1.6. The subsample of objects fitted by a star forming galaxy template 
has an even harder stacked spectrum, with V ~ 1.2 — 1.3. The obscured fraction is larger at lower fluxes, lower redshifts 
and lower luminosities. X-ray absorption is less common than "optical" obscuration and its incidence is nearly constant 
with redshift and luminosity. This implies that at high luminosities X-ray absorption is not necessarily related to optical 
obscuration. The estimated surface densities of obscured, unobscured AGN and type 2 QSOs are respectively 138, 59 
and 35 deg" 2 at F > 10" 14 erg cm" 2 s" 1 . 

Key words. X-ray - AGN - X-ray surveys 



1. Introduction 

One of the goals of deep X-ray surveys is to probe the origin 
of the X-ray background (XRB) to the faintest flux levels; 
they have resolved into discrete sources more than 90% of 
the XRB in the 0.5 - 2 keV band and up to 80 - 90% in the 
2 — 10 keV band (iGiacconi et aTl l2002t lAlexander et all . 
l2003HMoretti et all 120031 : iDe Luca fc Molendil . l2004D . The 
resolved fraction of the XRB drops however to ~ 60% above 
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- 6 keV and - 50% above - 8 keV (|Worslev et all I2005D . 
Most of the sources detected in deep, pencil-beam X-ray 
surveys are characterized by poor counting statistics, pre- 
venting a detailed analysis of X-ray spectral properties, and 
optical counterparts are often too faint for spectroscopic 
follow-up. Medium deep surveys, covering larger areas, are 
useful to bridge the gap between known X-ray source pop- 
ulations at low redshifts and those required to model the 
background and to collect a large number of sources for 
which X-ray and optical spectral analysis are fe asible (e.g. 
iPiconcelli etail 120031: iGeorgakakis et all l2006fh They are 
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also more effective th an deep surveys to find rare objects, 
such as type 2 QSOs (|Fiore et all . l2003fh 

The XMM Medi um Deep Survey (XMDS, see 
IChiappetti et all 120051 hereafter Paper I) consists of 19 
X-ray pointings, of nominal exposure of 20 ksec, cov- 
ering a contiguous area of about 2.6 deg 2 . It also lies 
at the heart of the lar ger, shallowe r XMM Large Scale 
Structure (LSS) Survey ([Pierre et all . 12004 iPacaud et all . 
I2006h . which will cove r ~ 10 deg 2 and i s principally de- 
voted to clusters study (jPierre et alll2006T ). Several surveys 
at different wavelengths are associat ed to the XMDS: the 
VIRMOS VLT Deep Survey fVVDS. lLe Fevre et al.l . l200l 
and the Canada - France - Hawaii Telescope Legacy Survey 
(CFHTLSfl in the op tical, the UKIRT Infrared Deep 
Sky Survey (UKIDSS, iDve et all 120061: lLawrence et ail 
I2006D in the near-IR and the Spitzer Wid e- Area InfraRed 
Extragalactic Legacy Survey (SWIRE, lLonsdale et all 
20030 in the mi d-IR. Radio observ ations performed at 
VLA at 1.4 GHz dBondi et all l2003h and at 325 and 74 
MHz ([Cohen et all |2003| ) alsocover the XMDS area. 

In Paper I the catalogue of XMDS sources detected 
(S/N > 4er) in at least one of five energy bands 0.3 — 0.5, 
0.5 - 2, 2 - 4.5, 4.5 - 10 and 2-10 keV within the VVDS 
field (area ~ 1 deg 2 ) was presented together with tentative 
optical identifications. The logN-logS distributions were de- 
rived for X-ray sources in the full XM DS area separately in 
the two bands 0.5 - 2 and 2 - 10 keV. iGandhi et all {2006) 
computed the logN-logS in the same energy bands for X- 
ray sources in the whole XMM - LSS area, finding results 
in agreement with those of Paper I. 

Here we consider a sample of X-ray sources selected in 
the "hard" ,2 — 10 keV, band, in order to investigate the 
populations of obscured and unobscured AGN and discuss 
their multiwavelength properties in a way unbiased by the 
intensity of the sources in soft X-rays. In order to take 
advantage of the best multifrequency coverage available, 
we consider sources in the VVDS area. Using the Spitzer 
data we construct the mid-IR/optical to X-ray spectral en- 
ergy distributions (SEDs) of the sources to estimate red- 
shifts and classify AGN into different categories according 
to the best fitting template. We use the term "optically 
unobscured" AGN for objects fitted by a type 1 AGN tem- 
plate, and "optically obscured" AGN for objects fitted by 
type 2 AGN or a star forming galaxy template, indicating 
that the optical-UV emission from the AGN is at least par - 
tially hidden. In a companion paper (|Polletta et all |2007| ) 
the templates used and the observed SEDs are presented 
in detail and their dependence on luminosity and absorp- 
tion is discussed. A comparison between X-ray and optical 
properties of AGN with optical spectroscopy in the whole 
XMM-LSS is in progress (Garcet et al., in prep.). 

Independently of the SED classification, the X-ray spec- 
tra and/or hardness ratios allow us to estimate absorption 
in the X-ray band, assoc iated to interven ing gas. In uni- 
fied models for AGN (e.g. lAntonuccil [l993f ) . obscuration by 
dust and absorption by gas are thought to occur in a dusty 
torus surrounding the AGN. The increasing complexity of 
properties shown by individual AGN has lead to a revision 
of this simple scheme proposing different regions around 
the AGN as sit e s of absorption at different wavelen gths 
(e.g. lElvisl [20001 : iKrongold et all [20071 lElitzurl l2006h . We 
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Fig. 1. Layout of the 19 XMM - Newton pointings of the 
XMDS (circles; the dashed circle marks field G12, which 
was not analyzed because of high background level, see 
Paper I) and of the associated surveys: VVDS (solid rect- 
angle), CFHTLS Deep (Dl, short dashed rectangle), and 
Wide (Wl, area on the right side of the vertical solid line), 
UKIDSS (long dashed rectangle) and SWIRE (area on the 
right side of the diagonal solid line). 



will therefore distinguish optically obscured and X-ray ab- 
sorbed AGN and examine separately their dependence on 
X-ray flux, redshift and luminosity. 

The paper is organized as follows: the multiwavelength 
data set is presented in Section [21 Optical and IR identifi- 
cations are discussed in Section [3] while the X-ray, optical 
and IR properties of the sample are derived in Section [4j 
The template SEDs and fitting process leading to the esti- 
mate of photometric redshifts and to the AGN classification 
are described in Section [5] The X-ray spectral analysis is 
presented in Section [5] X-ray spectra are analyzed individ- 
ually for sufficiently bright sources, while for faint sources 
absorption is estimated from hardness ratios. A stacking 
technique is used to derive average X-ray spectra of sub- 
samples of AGN with different SED classification. The sur- 
face density of optically obscured and unobscured AGN and 
of type 2 QSOs is derived in Section [71 Section [5] is devoted 
to the comparison of the fractions of optically obscured or 
X-ray absorbed AGN as a function of redshift, X-ray flux 
and luminosity. Finally, in Section [5] the results are sum- 
marized. 

Throughout the paper Ho = 70 Km s _1 Mpc -1 , Ha 
0.7 and Ojvf = 0.3 are assumed. 



2. The multiwavelength data set 

The layout of the XMDS observations is shown in Fig. [TJ 
superposed are the areas covered by the various associated 
surveys at other wavelengths, VVDS and CFHTLS in the 
optical, UKIDSS in the near-infrared and SWIRE in the 
mid-infrared. Shallower XMM - Newton pointings in the 
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Fig. 2. 2 - 10 keV flux distribution for the 136 X-ray sources 
in the 3 a hard sample. 



context of the XMM - LSS lie all around the XMDS area 
(see Pierre et al., in prep.) 



2.1. X-ray data 

The sample includes 136 X-ray sources detected at > 3cr in 
the 2 — 10 keV band (the 3 o hard sample hereafter), which 
fall within the sky area covered by the VVDS photometric 
survey (~ 1 deg 2 , solid rectangle in Fig.[I}. This area ben- 
efits from the richest multiwavelength coverage as evident 
from Fig. [TJ 

The sources were extracted from the XMDS catalog de- 
scribed in Paper I, to which we refer also for details on 
the XMM - Newton observations and data reduction. For 
all the sources count rates and fluxes were obtained inde- 
pendently in 5 energy bands: 0.3 — 0.5, 0.5 — 2, 2 — 4.5, 
4.5 — 10 and 2—10 keV. Fluxes were computed for a sim- 
ple power law spectrum with spectral index T = 1.7 and 
the average galactic colu mn density in the XMDS region 
(Nff = 2.61 x 10 20 cm~ 2 . iDickev fc Lockmanl 1 1990( 1 sepa- 
rately in each energy band. Fig. [2]shows the 2 — 10 keV flux 

distribution; the lowest flux that we sample is ~ 10 -14 erg 

— ? — i 
cm s . 

Hardness ratios were computed for all sources (see sub- 
section [IT!]). For 55 sources we detect a sufficient number of 
net counts (> 50 in the 2 — 10 keV band) to attempt a spec- 
tral analysis for each source (see details in subsection 16. 1| . 



2.2. Optical and near-infrared data 

Broad band BVRI photometr ic observations from the 
VVDS (|McCracken et all 120031 ) are available for XMDS 
sources over an area of about 1 deg 2 . This photometry was 
obtained at the Canada France Hawaii Telescope (CFHT) 
with the CFH12K camera at limiting magnitudes of Bab ~ 
26.5, V AB ~ 26.2, R AB ~ 25.9 and I AB ~ 25.0 (50% com- 
pleteness for point sources). 



The U band imaging was performed over an effective 
area of ~ 0.71 deg 2 with the Wide Field Imaging (WFI) 
mosaic camera on the ESO MPI 2.2 m telescope at La Silla, 
Chile. Two different U filters were used, the ESO Z7/360 
filter and the Loiano Observ atory U filter. The li miting 
magnitude is Uab ~ 25.4 ( see iRadovich et ail . l2004h . 

A small area (~ 165 arcmin 2 ) within the VVDS was also 
observed in the J and K bands down to a limiting magni- 
tude of Jab ~ 24.2 and Kab ~ 23.9 (50% completeness for 
point sources) with the SOFI instrument mounted on the 
ESO NTT telescope. A detail ed description of th e K band 
imaging survey is reported in llovino et al.l (|2005l l . 

Optical spectroscopy with the Visible Multi Object 
Spectrograph (VIMOS) on the ESO - VLT UT3 in the 
VVDS area is still in progress; the project aims at observ- 
ing a representative large subsample of objects down to a 
limiting mag of Iab < 24 (jLe Fevre et alll2005l) . 9 sources 
in the present sample have been observed and for 8 of them 
a redshift has been derived. Additional spectroscopic red- 
shifts for 22 sources have recently been obtained from 2dF 
and VLT FORS2 observations performed in the context of 
XMM-LSS follow up campaigns (see Garcet et al., in prep.). 
Three other redshifts are from the literature. 

The XMDS area lies within the sky region covered by 
the CFHTLS, a large collaborative project between the 
Canadian and French communities. Observations use the 
wide field imager MegaPrime equipped with MegaCam, in 
the u*g'r'i'z' filters. Both the "Wide" survey field Wl (8 
deg x 9 deg) and the "Deep" survey field Dl (1 deg x 
1 deg) cover th e XM DS region at a depth of i' = 24.5 
(iHoekstra et all |2006[) and i' = 26.1 (50% completeness 
limit. ISemboloni et all 120061) . respectively. In the following 
we will use the Dl notation for data from the CFHTLS 
Deep and the Wl notation for data from the CFHTLS 
Wide. Wl observations have been recently completed; the 
coverage shown in Fig [1] is deduced from data available to 
us at the time of analysis. 

Near-infrared observations of the XM DS area are also 
m progress m the cont ext of the UKIDSS (|Dve et all 120061 : 
lLawrence et~aT1 . I2006D . The survey uses the Wide Field 
Camera (WFCAM) of the 3.8 m United Kingdom Infrared 
Telescope (UKIRT). The XMM-LSS region is one of the 
four target fields of the Deep Extragalactic Survey (DXS). 
Observations in the J and K filters down to J = 22.3 and 
K = 20.8 (Vega system) are in progress. About 0.8 deg 2 of 
sky have been observed up to now and part of t he data are 
availa ble in the UKIDSS Early Data Release (jDve et all 
2006). Th e photometric system used in the UKIDSS is de- 
scribed in lHewett etail (|2006( l. 



2.3. Mid-infrared data 

The XMDS an d XMM-LSS fields are covered by the 
SWIRE survey (|Lonsdale et all l2003f ). Observations were 
performed with the Infrared Array Camera (IRAC) at 
3.6, 4.5, 5.8 and 8.0 /j,m and with the Multiband Imaging 
Photometer (MIPS) at 24, 70 and 160 ^m to a 5ct depth of 
4.3, 8.3, 58.5, 65.7 ^Jy and 0.24, 15, and 90 mJy, respec- 
tively. The area covered in the XMDS region is about 2.5 
deg 2 , corresponding to about 95% of the field. 12 sources 
in our selected sample are outside the region covered by 
SWIRE. De tails on the SWIRE data and source catalogs 
are given in ISurace et all (|2005l ). 
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2.4. Radio data 

There are also radio observations associated with the 
XMDS: the VLA VIRMOS Survey, which covers the VVDS 
area at a de pth of 80 /xJy (5a limit) and a re s olutio n of 6" 
at 1.4 GHz (|Bondi et all 120031 : ICiliegi et all . I2005D and a 
low frequency radio survey performed for the XMM-LSS 
also at VLA, which covers 5.6 deg 2 at a depth of 4 mJy at 
74 MHz and 110 de g 2 at a depth of 275 mJy at 325 MHz 
(jCohen et al.l . l2003f) . 

3. Optical and Infrared identifications 

Most (80%) of the 136 X-ray sources in the present sam- 
ple were already included in the 4cr catalogue presented in 
Paper I since they are also detected at > 4er in the softer 
bands. 24 X-ray sources are considered here for the first 
time. Although most sources had already been assigned op- 
tical counterparts, we have repeated the identification pro- 
cedure on the whole sample in a semi automatic way, that 
takes into account the experience accumulated in Paper I 
and the CFHTLS and SWIRE data now available. 

Access to the whole VVDS and CFHTLS catalogues 
and images is restricted: photometric data and positions are 
provided only for optical sources within a fixed radius from 
the X-ray positions. The same is true for the SWIRE data. 
We then associated to each X-ray source all combinations of 
optical and IR objects in the considered catalogues within 
a search radius of 6". Objects in the VVDS, CFHTLS and 
SWIRE catalogues were matched only a posteriori. 

We computed the probability of chance coincidence 
between an X-ray source and all optical VVDS, optical 
CFHTLS and infrared SWIRE candid ates within the search 
radius using the following equation (|Downes et all Il986t 
see also Paper I) 

p = 1 — exp(—n n(< m) r 2 ) 

where r is the distance between the X-ray source and the 
proposed counterpart (with a lower value fixed at 2", which 
roughly corresponds to the XMM - Newton astrometric un- 
certainty), and n(< m) is the density of objects brighter 
than the magnitude m of the candidate counterpart. We 
used / magnitudes for VVDS sources and i' magnitudes 
for CFHTLS objects (Dl data where available). We used 
the density n(> F 3 , 6 ) for infrared candidate counterparts. 
For each candidate counterpart there are therefore from 1 
to 3 values of p, depending if the object is detected in the 
VVDS, CFHTLS and SWIRE. We classified the probabili- 
ties as "good" (p < 0.01), "fair" (0.01 < p < 0.03) or oth- 
erwise "bad" and took as identification the object with the 
best probability combination. All tentative identifications 
were then checked by visual inspection using the VVDS 
finding charts. As reported in Paper I, astrometrical cor- 
rections were already applied to the X-ray fields, so we find 
again that 98% of the counterparts are within 4" from the 
X-ray corrected position, justifying our conservative choice 
of 6" radius. 

The probability criterion allows us to prefer one can- 
didate counterpart in the majority of cases, giving us 126 
secure identifications (out of 136 sources). Of these, 3 have 
only IR counterparts (i.e. no optical counterparts are de- 
tected down to Rab = 25.3) and will be referred to as 
optically blank fields. We notice that all sources covered by 
SWIRE (i.e. all but 12) are also detected in the IR. 



For the remaining 10 sources, the identification process 
is ambiguous leaving two or more possible counterparts, 
with similar probabilities. However, in 6 cases, the coun- 
terparts have similar magnitudes and colors allowing us to 
include these sources in parts of the discussion not involv- 
ing the redshift determination or SED classification. In the 
other 4 cases, the sources are completely dismissed. 

The X-ray, optical and infrared properties of sources of 
the 3 a hard sample are reported in Table [A~T1 For brevity, 
not all data used in this work are reported in the Table. 
The SWIRE catalogue is available through IRSA/Gator 
(http : / / irsa. ipac . caltech. edu/applicat ions/Gat or). 
We plan to publish the Catalogue of all XMDS X-ray 
sources with optical and IR identifications in a future 
paper. 

We also searched for UKIDSS counterparts of X-ray 
sources using a radius of 4", finding a near infrared coun- 
terpart for 72 X-ray sources. Generally UKIDSS sources 
are coincident with optical counterparts. There are however 
two exceptions: source XMDS 449J for which the UKIDSS 
source lies between the two possible counterparts, at a dis- 
tance of about 3" from both, and XMDS 760, for which 
there are two possible UKIDSS counterparts, both within 
~ 1" from the optical counterpart. The first case is one 
of the 4 X-ray sources that we could not identify, and the 
UKIDSS detection did not allow us to resolve the ambigu- 
ity. In the second case we associated to the optical counter- 
part the brightest UKIDSS source. 

33 X-ray sources in the 3cr hard sample (24%) have a 
radio counterpart at 1.4 GHz, one of them is also detected 
at 325 MHz. One is however associated with the spectro- 
scopically confirmed cluster XLSSC 025. We will not use 
the radio information in this work, but we checked the con- 
sistency of the radio fluxes with the templates used to fit 
the optical and infrared spectral energy distributions of our 
objects (see Section [5]). The correlat ion between th e X-ra y 
and radio luminosities is explored in lPolletta et al.1 (|2007l ). 



4. The AGN sample 

Ignoring two sources corresponding to spectroscopically 
confirmed galaxy clusters (XLSSC 025 and 041, see 
iPierre et alll200"a ). we computed the X-ray to optical ra- 
tio for 124 sources with secure identifications, using the 
equation given in section 6.2 of Paper I. For about 20 ob- 
jects for which VVDS magnitudes were unreliable because 
of saturation or unfavorable position in the field of view, the 
CFHTLS r' band magnitudes wer e used, with the appro- 
priate conversion factor taken from lSilverman et~a l. (2005). 
The X-ray to optical ratio is shown as a function of X-ray 
flux in Fig. El 

About 80% of the sources fall in the typical range 
of X-ray to optical ratio correspo nding to the locus of 
AGN (0.1 < Fx/Fr < 10 , see e.g. lAkivama et "all . 120001 : 
iHornschemeier et al.l 120 01). while about 20% of the sources 
have Fx/Ffj > 10, which correspond s to heavy absorption 
in th e optical and/or high redshift (jHornschemeier et all 
I2001h . This issue will be further developed in subsection l7.ll 



3 For brevity, in the text we label single sources with their 
XMDS identification number. The names of the sources, com- 
plying the IAU standard, along with the associated identifiers, 
are reported in Table \AA\ 
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Fig. 3. X-ray to optical ratio (2—10 keV band vs R band) as 
a function of X-ray flux for sources in the 3 a hard sample. 
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Fig. 4. B — I distribution for optical counterparts of X-ray 
sources in the 3 er hard sample. The dotted line marks the 



Diagonal dotted lines indicate loci of constant R magnitude division between the adopted definition of "blue" and "red" 
while horizontal solid lines mark the region of canonical objects (see text). 
AGN (0.1 < F x /F R < 10). Lower limits mark the optically 
blank fields, whose R magnitu de was fixed to 25.3 (see Fig. 
13 in lMcCracken et all . 120031) . 15 



Only two sources fall significantly below the AGN bor- 
derline (XMDS 1248 and 842): both appear extended in the 
optical as well as in the infrared images. The first (XMDS 
1248) has a low hardness ratio, consistent with no intrinsic 
absorption in the X-ray spectrum, while the second (XMDS 
842) has a higher hardness ratio, possibly indicating X-ray 
absorption (N# ~ 10 22 cm" 2 ). Using photometric redshifts 
(see Section [S]), we obtained X-ray luminosities of ~ 10 41 
erg s _1 for both of them, even after correcting for absorp- 
tion (see Table |A~T|) . We classify both provisionally as nor- 
mal galaxies, though we can not exclude the presence of 
a low luminosity AGN or even a Compt on thick AGN in 
XMD S 842 (see e.g. FSC 1021+4724 in lAlexander etafl 
2005). Another source in the sample has Lo.5-10 < 10 
erg s _1 , XMDS 178, however its X-ray to optical ratio of 
0.27 is in the typical AGN range. On the basis of the X-ray 
to optical rati o , we r etain this source in the AGN class. In 
iPolletta et alj (|2007f) slightly different criteria are adopted 
for these borderline objects. 

To summarize, on the basis of the X-ray to optical flux 
ratios 122 X-ray sources with unambiguous identification 
can be classified as AGN and 2 as normal galaxies. In the 
following subsections we will discuss the optical and IR 
properties of this sample. 



4.1. Optical magnitude and colors 

The B — I color and the Rab magnitude distributions for 
the identified sources are shown in Fig. [4] and [5l respec- 
tively. 

The B — I color distribution shows a high peak at 
B — I < 1.0, and a tail extending up to B — I ~ 4. Based 
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Fig. 5. R magnitude distribution for the 3 a hard sam- 
ple. Solid histogram: total sample; dotted histogram: blue 
sources (see text); dashed histogram: red sources (see text). 



on the observed color distribution we adopt a somewhat 
arbitrary threshold of B — I = 1.0 to divide the sample 
into two roughly equal size samples of "blue" objects, with 
B — I < 1.0 (43% of all sources), and "red" objects, with 
B — I > 1.0 (57%). As will be shown later, this criterion, 
although crude, proved to be a good one for a rough sep- 
aration between type 1 (i.e. broad line) AGN and type 2 
(narrow line) or star forming galaxy-like AGN based on 
observed quantities alone, and is substantially confirmed 
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Fig. 6. B — I color vs I magnitude for sources in the 3 a 
hard sample. Empty circles are blue objects, filled circles 
are red objects. Overimposed to points are the evolutionary 
tracks for an Sc galaxy (solid line), a type 1 QSO (dotted 
line) , a type 2 AGN (short-dashed line) and a type 2 AGN 
plus a starburst component (long-dashed line). The B band 
absolute magnitudes assumed are —15.7 for the Sc galaxy 
template, -22.3 for the QSOl, -16.9 for the type 2 AGN 
and —23.9 for the type AGN plus starburst. 



Fig. 8. IRAC color - color plot for sources in the 3 a hard 
sample. Empty circles are blue sources, filled circles are 
red sources, crosses are the optically blank fields or sources 
with undefined color classification. Left pointing arrows are 
sources undetected in the 5.8 /xm band, down pointing ar- 
rows are sources undetected in the 8.0 /im band. 5 a upper 
limits are used. The d ashed lines mark t he region expected 
for AGN according to lLacv et all (|2004l ) 



by the more detailed (but model-dependent) classification 
based on the spectral energy distributions. 

The magnitude distributions of these two broad classes 
are plotted separately in Fig. [5] on average, blue sources 
(B-I < 1.0) are brighter, with a peak at R ~ 20 and 90% 
of objects at R < 22, while red sources (B — I > 1.0) have 
a broader distribution, extending from R ~ 18 to R ~ 26. 

In Fig. [5] we show the B — I color as a function of 
the I magnitude for our sources, along with the evolu- 
tionary tracks for various templates: a late spiral galaxy 
(Sc, solid line), a type 1 QSO (dotted line), a type 2 AGN 
(short-dashed line) and a type 2 AGN pl us a starburst 
comp onent (long-dashed line; see below and lPolletta et all 
l2007f h The effects of absorption due to the Intergalactic 
Medium (IGM) have been tak en into accoun t at high red- 
shift (z > 2.5) as prescribed in lMadaul jl995). Blue sources 
are near the QSOl track, while red objects are generally 
consistent with star forming galaxies and AGN2 tracks. 
However for magnitudes fainter than Iab — 23 the dif- 
ferent track cross, and type 1 AGN, type 2 AGN and star 
forming galaxies have similar colors. 

4.2. X-ray to infrared ratios 

In Fig. [Jj we plot the ratios of X-ray to infrared fluxes at 
3.6/im (left panel) and 24/im respectively (right panel) as a 
function of the X-ray flux. Sources are all clustered in the 
same region with no clear separation between blue and red 
sources. 



Two typical loci of local sources are shown in the right 
panel of Fig. \7\ the area at -0.8 < log(F x /F IR ) < is 
the region occupied b y hard X-ray selected AGN (from 
IPiccinotti etaTl 11982ft with IR emission and z < 0.12; the 
area close to \og{Fx/ Fm) = — 3 is the reg i on oc cupied by 
local starburst galaxies fro m Ranalli et al.l (|2003f ). adapted 
from lAlonso-Herrero et ail (|2004l ). No objects with X-ray to 
infrared ratios typical of local starburst galaxies are found 
in our sample. 80% of the objects have X-ray/infrared ra- 
tios — 1 < log(Fx/i/ 3 . 6 S , 3.6) < (i.e. within a factor of 10), 
and 98% of them have -1.2 < \og{F x / v 3i6 S M ) < 0.6 (i.e. 
a factor of ~ 70). The most discrepant object is one of the 
two normal galaxies with low X-ray to optical ratio (see 
above). The X-ray to optical ratios for the same sources 
ranges from ~ 0.1 to ~ 60 (i.e. a factor of 600 excluding 
lower limits, see Fig. This implies that the IR flux is a 
better diagnostics of the X-ray flux compared to the opti- 
cal, a behaviour likely due to the smaller extinction in the 
IR and to the fact that nuclear light absorbed by dust is 
likely re-radiated in the IR. 

The observed range in the Fx/v-zaSza plot is fully 
cons istent with other X-ray and 24 urn samples , 
fe.g. lAlonso-Herrero et all 120041 : iFranceschini et all 120051 : 
iPolletta et al.l . l2006f ), but broader than that of local hard X- 
ray selected AGN of IPiccinotti et all (| 19821 ). This broader 
dispersion is not surprising given th e better sensitivity of X- 
ray observations with respect to the IPiccinotti et all £l982) 
data. 

A broad range in the X-ray to infrared ratio could 
be caused by different amounts of absorption in different 
sources that depresses the observed X-ray flux, but not the 
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Fig. 7. X-ray to IR ratio as a function of X-ray flux for sources in the 3<r hard sample. Empty circles are blue sources, 
filled circles are red sources, crosses are the optically blank fields or sources with undefined color classification (because 
B or / magnitudes are not available). In the left panel IRAC 3.6 fim flux is used, while in the right panel MIPS 24 fim 
is used. Lower limits are sources detected in one or more IRAC bands and undetected in the MIPS 24 /jm band, where 5 
a upper limit is u sed. The diagonally sha ded area is the region occupied by hard X-ray selected AGN with IR emission 
and z < 0.12 from l Piccinotti et al.l (|1982f ); the horizontally shaded area is the region occupied by local starburst galaxies 
from lRanalli et al.1 (|2003| ). 



infrared emission. Alternatively, it could be an intrinsic dis- 
persion in the AGN SEDs that is not sampled properly in 
local objects. If this dispersion were due only to absorp- 
tion in the X-rays, it would imply a broad range of column 
densities, up to 1.5 x 10 23 cm -2 , consistent with the distri- 
bution of measured column densities (see subsection 16. 2|l . 
However, the similarity in the distribution of flux ratios of 
blue and red sources is not observed in the column density 
distribution, the majority of blue sources being unabsorbed 
and the majority of red sources being absorbed. These argu- 
ments suggest that the variety of the intrinsic SED shapes 
that characterize the AGN population is a more likely ex- 
planation and that such a variety is also observed for opti- 
cally blue AGN. I n fact a recen t study of X-ray and 24^tm- 
selected AGN by iRigbv et ail (|2005f ) shows that there is 
no correlation between the ratio F(X) / ^24-F(24^m) and the 
amou nt of absorp t ion in the X-rays, or their optical prop- 
erties. |Mvis^r^L| (|1994| ) measure a dispersion of a factor of 
10 at 24/im for a large sample of optically-selected quasars 
after normalizing their SEDs at 1 /im, consistent with the 
observed dispersion in the X-ray/infrared flux ratios of our 
sample. An analysis of the SEDs of the AGN in the sam- 
ple is presented in t he next Section and in more detail in 
IPolletta et ail (|2007[ ). 

IRAC infrared colors proved to be a useful di agnostics to 
identi fy AGN among IR sources; in particular, lLacv et all 
(|2004h found that the 8.0/4.5 fim ratio vs the 5.8/3.6 /zm 
ratio plot is effective in isolating AGN in IR selected sam- 
ples, which have red colors (i.e. high values of the ratios) 
in bot h axis. In Fig. [8] we reproduce the plot of lLacv et al.l 
(2004) for sources in our sample, and we find that the vast 
majority of them (both optically blue and red) lies in the 
region expected for AGN. At the boundaries of the AGN re- 
gion there could b e contamination by low redshift galaxies 
(|Lacv et al.l . [2004) ; in fact, all the objects near the borders 
of the AGN region in Fig. [8] have a red optical color. The 
AGN with the reddest IR colors are predominantly blue in 



the optical, while optically red AGN show a broad range of 
IR colors. 



5. Photometric redshifts and SED classification 

Taking advantage of the excellent multiwavelength cover- 
age from the optical (VVDS, CFHTLS), to near- and mid- 
infrared (UKIDSS and SWIRE) we constructed broad band 
SEDs for all the 124 identified sources. We then fitted the 
observed SEDs (taking into account also upper limits) with 
various templates in order to determine photometric red- 
shifts. We used 20 templates that represent normal galax- 
ies (11: 1 elliptical, 7 spirals and 3 starbursts), composite 
galaxy + AGN (3: starburst + AGN) and AGN (6: 3 type 1 
AGN, 3 type 2 AGN) and cover the wavelength range from 
1000 A to 500 fim. These were derived from the observed 
SEDs of objects representing the different classes. The ef- 
fects of dust extinction were taken into account by red- 
dening the reference templates according to the extinction 
curve derived in high redshift starbursts by ICalzetti et al.l 
()2000D . In order to limit degeneracies in the best fit solu- 
tions we limited the extinction Ay to be less than 0.55 
mag and included templates of highly extincted objects 
to fit more heavi l y obs cured sources. The HYPERZ code 
(jBolzonella et al.1 . |2000| ) was used to fit the SEDs and find 
the best-fit solution. A full description of the templates and 
a detailed discussion of the SED fitting proc edure and pho- 
tomet ric redshift estimates are presented in IPolletta et al.l 
((20071) . 

A number of spectroscopic redshifts are available to as- 
sess the quality of our photometric redshift determination. 
For 22 objects redshifts were obtained in the context of 
the XMM - LSS follow up programs and made available 
to us (Garcet et al., in prep . ). Red shifts for two sources 
were taken from lLacv et ail (|2006), who present optical 
spectroscopy of luminous AGN selected in the mid-IR from 
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Fig. 9. Photometric vs spectroscopic redshifts. Solid line 
is the z p h = z sp relationship, dotted lines mark the 20% 
"error" in 1 + z p h and dashed lines mark the 30% "error" 
in 1 + Zph- Hexagons are type 1 AGN, empty squares are 
type 2 AGN and asterisks are SGFs, as characterized by 
their SEDs discussed in subsection l5.ll 



Spitzer observations. For XMDS 842 a redshift is available 
from NEE0. 

The VVDS spectroscopic sample (see Gavignaud et al. 
2006 for type 1 AGN) yields redshifts for 8 sources in the 
present sample. To obtain a larger redshift comparison set, 
we added 16 additional sources from the larger X-ray sam- 
ple discussed in Paper I having a redshift from the VVDS 
spectroscopic survey. For the latter similar photometric 
data are available so that photometric redshifts could be es- 
timated with the same procedure described above. In total, 
the spectroscopic comparison sample consists of 49 sources. 
For 3 of them, falling outside the area covered by SWIRE, 
only optical data were available for the SED. 

Photometric and spectroscopic redshifts are compared 
in Fig. [5J The reliability and accuracy of the photomet- 
ric redshifts are usually measured via the fractional error 
Az = ^ z p^~ z "p" ^j an( } the rate of catastrophic outliers, 

defined as the fraction of sources with \Az\ > 0.2. For our 
49 objects, the mean Az is consistent with 0.00, with al a 
dispersion of 0.12, and the outlier fraction is 10%. These 
results are significantly better than previously obtained 
for AGN samples, where the fraction of outliers is usually 
highe r than 25% (|Kitsionas et all 120051 : iBabbedge et all 
120041) . The achieved accuracy still does not allow us to 
consider photometric redshifts as fully reliable for individ- 
ual sources, however it is adequate for a statistical analy- 
sis of the population . For a more detailed discussion, see 
iPolletta et all (|2007t ). 

The distribution of the 124 photometric redshifts (in- 
cluding the optically blank fields, for which only IR fluxes 
were used) is shown in Fig. [10] The majority (60%) of 

4 http : //nedwww . ipac . caltech . edu/ 



sources has z < 1, with a tail extending up to z ~ 4. These 
results are consistent with the redshift distribution of other 
X-ray selected samples with similar or deeper flux limit (e.g. 
IBarger et all l2003t iHasingerL 120031 : IBarger et all 120051 ). 

5.1. Spectral energy distributions and classification 

According to the template which gives the best-fit solution, 
we assigned sources to one of the following broad classes: 
type 1 AGN, type 2 AGN, or SFG. The type 1 AGN class 
corresponds to sources best-fitted with a QSOl template. 
The type 2 AGN class includes sources best-fitted with ei- 
ther the Seyfert 2 templates, or the composite AGN + star- 
burst templates, or the QS02 template. The SFG class in- 
cludes sources fitted by a spiral or a starburst template. 
Elliptical templates never yielded best fit solutions. 

Examples of observed SEDs with their best fit templates 
are presented in Fig. llH For sources with both optical, near 
and mid-IR data, the photometric classification should be 
reliable since the SED shape of the different classes has spe- 
cific signatures that can be easily identified. Interestingly, 
while photometric redshifts for type 1 AGN might be the 
most uncertain, their classification is instead rather easy. 
Note however that the Seyfert 1.8 template appears inter- 
mediate between type 1 and type 2 AGN (see Fig. [TTj) . 
There is a large variety of SED shapes among the tem- 
plates used for type 2 AGN, composite and star forming- 
galaxies. In case the fit is not optimal or when only few IR 
data points are available, the separation between the var- 
ious classes is uncertain as can be guessed comparing the 
SEDs in the left panel in Fig. [TTJ 

The SED fitting procedure yields 39 type 1 AGN (32%), 
61 type 2 AGN (49%) and 24 SFG (19%). 

Comparing the SED classification with the spectro- 
scopic one, we find that all the 16 objects classified as type 
1 AGN from the fitted template indeed show broad emis- 
sion lines in their optical spectra. Thus a photometric type 
1 AGN classification appears unambiguous. 

On the other hand, there are 10 objects spectroscopi- 
cally classified as type 1 AGN, which are instead not recog- 
nized as such by the SED fitting procedure, indicating that 
our method systematically underestimates the fraction of 
broad line AGN. Specifically of the 10 misclassified objects 
8 are fitted by a Seyfert 1.8 template (all with Ay close to 
0), one by a QSO 2 template and 1 by a SFG template. 
These objects appear to be dominated by star-light emis- 
sion in the optical and near-IR where the AGN continuum 
does not emerge clearly, although broad emission lines are 
visible in the optical spectrum. Of the remaining 23 objects 
without broad lines in their optical spectra only 5 are fitted 
with a Seyfert 1.8 template, in three cases with significant 
extinction. We conclude that SEDs fitted by Seyfert type 
1.8 templates are intermediate between type 1 and type 2 
objects and that our method systematically underestimates 
the objects spectroscopically classified as type 1. 

The sources photometrically classified as SFGs do not 
show any AGN signature at optical and IR wavelengths, 
however the X-ray to optical and X-ray to IR ratios and 
the X-ray luminosity unambiguously point to the presence 
of AGN activity also in these objects. 

In the following we will define optically "unobscured" 
AGN all sources fitted by a type 1 AGN template. These 
sources are expected to unambiguously correspond to broad 
line AGN. We will define all other sources (i.e. having ei- 
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Fig. 10. The distribution of photometric redshifts is shown as solid histogram in both panels. The dotted and dashed 
histograms (left and right panels, respectively) refer to the subsamples of unobscured and obscured AGN, respectively, 
discussed in subsection 15. II 
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Fig. 11. Observed SED (filled circles) and redshifted best-fit templates (grey solid curves) of 13 sources: 2 with star- 
forming like SEDs (left panel), 5 with type 2 AGN SEDs (left panel), and 6 with type 1 AGN SEDs (right panel). 
Downward pointing arrows correspond to 5rj upper limits. The source sequence number and best-fit template names are 
listed in the same order as the SEDs are plotted. 



ther type 2 AGN or SFG like SEDs) as optically "obscured" 
AGN. As shown above, the latter group may include some 
AGN with broad emission lines, but with a SED dominated 
by the host galaxy in the near-IR. We will take into ac- 
count where relevant that the number of unobscured ob- 
jects should be corrected upwards by a factor 1.6 (and the 
number of obscured objects reduced accordingly). 

In Fig. [T21 we compare the classification based on the 
SED shape with the optical color B — I as a function of red- 
shift. The horizontal dashed line corresponds to the thresh- 



old between blue and red sources (B — I = 1). 90% of the 
unobscured AGN have blue optical color and 84% of ob- 
scured AGN are red. Thus the simple classification based 
on observed color appears in retrospect rather successful 
when compared with the more sophisticated template fit- 
ting procedure. However, while at z < 1.6 the SED and 
"color" classifications practically coincide (except for 2 ob- 
scured objects near the borderline), at larger redshifts there 
is a degeneracy among the different evolutionary tracks so 
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Fig. 12. B — I distribution of sources in the 3 a sam- 
ple as a function of photometric redshift. Empty diamonds 
are unobscured AGN, filled diamonds are obscured AGN. 
Overimposed to points are the evolutionary tracks for an 
Sc galaxy (solid line), for a type 1 QSO (dotted line), for a 
type 2 AGN (short-dashed line) and a type 2 AGN plus a 
starburst component (long-dashed line). 

that the optical color alone is not indicative of a spectral 
type. 

Based on the photometric classification, the redshift dis- 
tributions of unobscured and obscured AGN can be derived. 
They are shown in the left and right panels of Fig. [TU1 
respectively. The two are clearly different, the first be- 
ing broader and reaching higher redshifts, while the sec- 
ond is more concentrated at z < 1 (72%). Several authors 
(e.g.lEckart et aUl2006t ISteffen et a ll 12004 Pfr eister et all 
2005; lLa Franca et all I2005Q find "different redshift distri- 
butions for type 1 and non type 1 AGN. An analysis of the 
fraction of unobscured and obscured AGN as a function of 
redshift will be discussed in Section [SI 



6. X-ray spectral properties 

We studied the X-ray spectral properties of our sample per- 
forming spectral fitting for individual sources with a suffi- 
cient number of counts fsubsection l6.1[) and a hardness ratio 
analysis for for faint sources to obtain individual values of 
Nh (subsection 16. 2[) . A stacking technique was then used 
to study systematic trends in the whole sample (subsec- 
tion 16. 3p . The general hardness ratio definition is 



HR 



CRh ~ CRs 



(1) 



CR H + CR S 

where CRh and CRs are the count rates in the hard and 
in the soft band, respectively. In subsection 16 . 21 we consider 
three hardness ratio values for each sources, HRcs, calcu- 
lated using the energy bands 0.5 — 2 (B) and 2 — 4.5 keV 

(C) , HRdc, using the energy bands 2 — 4.5 and 4.5 — 10 keV 

(D) and HR, computed between the energy bands 0.5 — 2 
and 2-10 keV. 



6.1. Individual sources 

We extracted X-ray spectra for all sources in the 3 a hard 
sample having at least 50 net counts in the MOS + pn 
merged image in the 2 — 10 keV band. There are 55 X-ray 
sources in our catalogue which satisfy this criterion; 24 of 
them are optically unobscured AGN, and 31 are optically 
obscured AGN. 

Counts were extracted for each source using the XMM 
- Newton Science Analysis System (SAS) evselect task 
in a circular region with a radius of 20", corresponding, 
for a point-like source, to an encircled energy fraction of 
~ 70 — 75% (off-axis angles between and 10'). The pn 
data were used, unless the source was close to a CCD gap, 
in which case we used the MOS data, fitting simultaneously 
MOS1 and MOS2. Background counts were extracted from 
the nearest source free region, excluding areas near gaps in 
the CCD array. We used the SAS rmfgen task to create 
response matrices (one for each camera and each XMDS 
pointing) and arfgen to generate ancillary response files 
(one for each source). 

X-ray spectra were analyzed using the XSPEC package 
(v. 11.3.1). We considered the energy range 0.3 — 10 keV. 
When the number of counts was large enough, data were 
binned in order to have at least 15 or 20 counts for each 
energy channel and y 2 statistic s was used, otherwise we 
used Cash statistics (|Cashl . [T979l) . which, however, does not 
give a "goodness of fit" evaluation, like the \ 2 . In order to 
better match the spectral resolution of the instruments, we 
binned the data of these sources with few counts using a 
fixed number of PHA channels before fitting using the Cash 
statistics. 

We first fitted the spectra using a simple power law 
model with galactic absorption com puted at the XMDS 
posit ion (Nfj = 2.6 x 10 20 cm -2 , iDickev &; Lockmanl 
1990), plus a component for intrinsic absorption at z = 
( XSPEC model : phab s*zphabs*pow with abundance table 
of lWilms etalll2000h . 

For all spectra for which \ 2 statistics can be used in 
the fit (22 sources), we set both intrinsic column density 
and photon index as free parameters. Spectral fit results 
are reported in Table IB. II in Appendix [B] The errors in 
Tables and Figures correspond to the 90% confidence level 
for one interesting parameter. The average photon index 
is r ~ 2.1 and N# < 21 cm ~ 2 . This is consistent with 
their location in the hardness ratio plot, where they cluster 
around BR C b = -0.5 and HR DC = -0.5 (see Table ED 
and subsection 16. 2|) . Therefore they cannot be considered 
representative of the whole sample. 

Since more than half of the X-ray spectra in the sam- 
ple do not have a sufficient number of counts to perform 
a fit with both T and N# free parameters, we fixed the 
photon index for all objects, in order to obtain an esti- 
mate of the column density. We used two different values 
of the photon index, T = 2.0 and T = 1.7, both appropri- 
ate fo r AGN (jTurner fc Poundsl Il989t iNandra k Pounds! . 
1994). Spectral fit results for the simple absorbed power 
law model for each source with both T = 2.0 and T = 1.7 
are reported in Table IB. 21 in Appendix [B] where we also 
list the sources with peculiar fits. The best fit values of N# 
obtained with r = 2.0 are higher than those obtained with 
photon index frozen to 1.7, by about Alog(N ff ) = 0.2. The 
two column density estimates are consistent within errors 
in 90% of cases. 
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Fig. 13. Intrinsic column density vs photometric or spec- 
troscopic redshift. Empty diamonds are optically unob- 
scured AGN, filled diamonds are optically obscured AGN. 



We will consider in the following only the distribution 
obtained fixing T = 2.0, except for two sources (XMDS 124 
and 779): in these two cases we were able to find a stable so- 
lution only fixing the photon index to T = 1.7. The column 
density distributions of optically obscured and unobscured 
AGN turn out to be different: 16% of unobscured AGN (3 
out of 19) have N# > 10 21 cm -2 , while more than 55% of 
obscured AGN (19 out of 34) have N H > 10 21 cm -2 . We re- 
call that these are lower limits to the column density values, 
since we did not yet introduce the redshift dependence. 

Finally, we introduced the photometric or spectroscopic 
redshift, when available, in order to compute the intrinsic 
column density. The photon index was left free when the 
X 2 statistics could be used, otherwise it was fixed to 2.0. 
Spectral fit results, along with redshifts, are reported in 
Table [531 in Appendix |B1 Again, the distributions for opti- 
cally obscured and unobscured AGN are different, with 63% 
of obscured AGN (19 out of 30) having N ff > 10 21 cm" 2 
and 36% (11 out of 30) with N H > 10 22 cm" 2 . For com- 
parison, only ~ 20% of unobscured AGN (4 out of 21) have 
N H > 10 21 cm" 2 and 10% (2 out of 21) have N H > 10 22 
cnr 2 . 

In Fig. [T3] the intrinsic column density is shown as a 
function of redshift; this figure is qualitativ ely consistent 
with those presented in other surveys (e.g. lEckart et all 
2006) and shows no obvious trend with z, although we also 
notice the paucity of high redshift sources with well con- 
strained measures of N# . 

6.2. X-ray absorption 

The two hardness ratios HRcb and ERdc defined above 
are compared in Fig. 1141 As expected, most sources lie along 
the values expected for a single power law model. We fur- 
ther distinguish obscured/unobscured sources with differ- 
ent symbols. Less than 10% of the objects classified as op- 
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Fig. 14. X-ray color - color plot for sources in the 3 a hard 
sample. Energy bands involved are: 0.5 — 2 (B), 2 — 4.5 (C) 
and 4.5 — 10 keV (D). Sources are classified as optically un- 
obscured AGN (empty diamonds) or obscured AGN (filled 
diamonds). We also mark the hardness ratios computed for 
a simple absorbed power law spectral model, with T = 2 
and logN# = 21, 22, 23 for z — (four point stars con- 
nected by the solid line, from left to right) and logN# = 
22, 23, 24 for z = 1 (five points stars connected by the 
dashed line, from left to right). 



tically unobscured AGN have HR CB > -0.3 (N H ~ 10 215 
at z — 0), while more than 40% of the sources classified as 
obscured AGN have HRcb > —0.3 indicating that X-ray 
absorption and an obscured classification are often associ- 
ated. 

For a quantitative estimate of the absorbing column 
N# we used the results of the spectral fits described in 
Section |6~T1 with T fixed to 2.0 for the 51 brightest sources 
and computed the column density from the X-ray hardness 
ratios in the remaining cases in the following way. 

We used the standard hardness ratio HR, computed be- 
tween the 2 — 10 and the 0.5 — 2 keV bands. We made 
simulations using XSPEC to obtain hardness ratios cor- 
responding to typical values of N# ranging from 10 20 to 
10 24 cm -2 . A simple power law model with photon index 
r = 2.0 was assumed, consistently with the model used 
for the X-ray spectral analysis. The simulations and the 
objects for which a spectrum could be extracted define a 
clear relationship between N# and HR for HR > —0.5, 
while, below these values, the HR - N# relation degener- 
ates. We therefore fixed the latter value, corresponding to 
Nh ~ 10 21 cm -2 , as a threshold below which all column 
densities are fixed to the galactic value. By interpolation we 
computed the observed Nh corresponding to the hardness 
ratio of each source. The intrinsic column density was then 
obtained from the observed one using the photometric (or 
spectroscopic, when available) redshift and the expression 
N^f r = N^ s (l + z) 2 6 (iBarger et alJ . l2002ft . also when the 
observed column density was estimated from the spectrum 
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Fig. 15. Distribution of observed (left panel) and intrinsic (right panel) column densities. Solid histogram refers to the 
whole sample, dotted and dashed histograms refer to optically unobscured AGN and optically obscured AGN, respectively. 



(i.e. in this cases we did not use the N^* r obtained by the 
XSPEC and reported in Table IB~2l but we recomputed it 
from the observed value, in order to be more consistent with 
the column density estimates obtained from HR) . W^ tr was 
not computed when the observed column density was fixed 
at the galactic value. 

The observed and intrinsic column density distributions 
are reported in Fig. [TSJ Different lines (dotted/dashed) re- 
fer to the Nh distribution for optically unobscured and ob- 
scured AGN respectively. Unfortunately, our choice of set- 
ting a fixed value for low Njj creates an artificial gap in 
the distribution. The majority of optically obscured AGN 
are X-ray absorbed (N^ l * r > 10 22 cm -2 ), as expected. 
We find that also 12 unobscured AGN (more than 30%) 
have Wfi tr > 10 22 cm -2 , ft is well known that N# values 
are less well constrained with increasin g redshifts ( see e.; 
lEckart et all . l2006t lAkvlas et all . l2006HTozzi et all 12006 
since the absorption cut - off shifts to lower energies and be- 
comes comparable to the galactic values or even drops out 
of the observed band. For z ^ 1.5, an intrinsic Nw of 10 22 
cm -2 corresponds to an observed column density 10 21 
cm -2 . Unobscured AGN, whose hardness ratios generally 
cluster around —0.5 (which corresponds to N^ s ~ 10 21 
cm -2 ), are more severely affected by this problem than ob- 
scured AGN, which have a broader HR distribution. We 
have partially compensated for this effect with our choice 
of fixing intrinsic columns to when the observed hardness 
ratio HR ~ —0.5. Moreover, we have at least 4 examples of 
optically unobscured AGN (XMDS 12, 258, 280 and 406) in 
which the observed N# is already larger than N# = 5 x 10 21 
cm -2 , ensuring that the large columns are not all to be at- 
tributed to the redshift effects. 

On the other hand, among the objects for which we are 
not able to estimate the column density (those with N# 
fixed at the galactic value), there could be some which could 
be really X-ray absorbed. In Fig. [16] we show the intrinsic 
column density of our objects as a function of redshift. The 
solid line shows the intrinsic N# values that would be de- 
rived at a given redshift, for an observed column density of 
10 21 cm -2 . Since the objects with HR < —0.5 should have 



N obs 



< 10 cm, , their intrinsic column density should lie 
below the solid line in Figure. It is therefore possible that 
we underestimate the number of X-ray absorbed objects 
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Fig. 16. Intrinsic column density vs photometric (or spec- 
troscopic, when available) redshift. Empty diamonds are 
optically unobscured AGN and filled diamonds are optically 
obscured AGN. The dashed line marks the threshold be- 
tween X-ray absorbed (Njj > 10 22 cm -2 ) and unabsorbed 
sources. The solid line shows the intrinsic column density 
that would be derived at a given redshift, for a measured 
column density of 10 21 cm -2 in the observer frame. 



for redshift z <; 1.5 (where the solid and dashed line cross). 
The column densities are therefore difficult to estimate at 
high redshift, but this should not affect our results. 



6.3. Stacking analysis 

Given that less than half of the sources in our sample have 
a sufficient number of counts to perform a spectral analysis, 
we used a stacking technique to measure the spectral prop- 
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Table 1. Mean spectral properties from the stacking anal- 
ysis of sources in the 3er hard sample detected in the pn at 
off-axis angle < 11'. N# is fixed to the galactic value. 



erties, averaged over the whole redshift range, of sources of 
different classification and in different flux intervals. 

For the stacking analysis we used only pn data, because 
of the pn larger effective area; we selected only sources 
which are not in or near a pn CCD gap or bad column. 
Moreover, since the pn point spread function (PSF) and 
the v ignetting a r e not well de t ermin ed for large off-axis an- 
gles (jGhizzardil . 120021 : iKirshl . l2006f ) we only used sources 
with off-axis angle 9 < 11'. This value allows us to obtain 
a significant number of sources (83), for which calibration 
should be still reliable. 30 sources are optically unobscured 
AGN and 53 are optically obscured AGN. 

We restricted the source area to a fixed radius of 15", 
which, on-axis, corresponds to 67% of the encircled energy 
for a point-like source. Since with this radius we sample the 
PSF core and not the wings, the encircled energy fraction 
has a weak dependence on the off- axis angle (at 9 = 10' 
67% of the encircled energy is within a 16" radius) and the 
energy dependence can also be neglected. Therefore, we can 
consider all sources together regardless of their position in 
the field. A background spectrum was extracted for each 
source in a circle of radius 80" in the nearest source free 
region. Ancillary response files were also produced for each 
source, while, as for single spectra analysis, we used one 
response matrix for each XMM - Newton pointing. 

The sample used in the stacking analysis covers a flux 
range from about 10 -14 to 1.2 x 10~ 13 erg cm -2 s _1 ; we 
divided it in 5 flux bins, chosen to have a sufficient number 
of counts in each bin (see Table [TJ on average, brighter 
sources have a greater number of counts, so in the higher 
flux bins a smaller number of sources is included). 

The spectra within the same flux bin were added using 
the mathpha task of f tools to produce a single spectral file. 
The same was done for background files. The auxiliary and 
response files were combined using the addarf and addrmf 
tasks of f tools, respectively. The combined spectra were 
grouped to a minimum of 20 counts per bin and were ana- 
lyzed using XSPEC. 

We fitted the stacked spectra in the 0.3 — 10 keV range 
using a single power law model with column density fixed 
to the galactic value. The fit results are reported in TableQ] 
The r values obtained for the whole sample are consistent 
within errors with T = 1.7 — 1.8 in all the flux bins. 

We then considered the optically unobscured and ob- 
scured AGN separately. We divided them in 4 bins, using 
two slightly different binnings for the two subsamples dic- 
tated by the available statistics. The results are reported 
in Table [5] and the photon index as a function of flux is 
shown in Fig. 1171 The difference between the optically un- 
obscured and obscured AGN populations is evident: for 
the unobscured AGN the measured photon index is con- 
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Fig. 17. Photon index obtained from the fit of stacked spec- 
tra as a function of X-ray flux for the optically unobscured 
AGN (empty diamonds), optically obscured AGN (filled di- 
amonds) and SFGs (crosses). Vertical bars are errors on T, 
while horizontal bars show the flux bin width. 



sistent with r = 2 over the whole flux range, while for 
the optically obscured it is consistent with T — 1.5 — 1.6. 
Therefore the observed average spectral slope of unobscured 
AGN is consistent with that of typical broad lin e AGN 
(jTurner k Poundsl 119891 : iNandra &; Pounds! . Il994[ ). while 
that of optically obscured AGN is significantly harder. No 
significant dependence of the spectral index with flux is 
fou nd for optically un o bscur ed or optically obscured AGN. 

iGeorgakakis et all (|2006l ) merged the X-ray spectra of 
hard X-ray sources detected by XMM - Newton at F 2 _s > 
2 x 10~ 14 erg cm~ 2 s _1 and having a n optical counter - 
part in the Sloan Digital Sky Survey (jYork et alj . 1 20001) 
with red color (g — r > 0.4). They found that the stacked 
spectrum of these sources has a spectral slope of T = 
1.47 ±0.04 (sources observed with t he THIN filter), consis- 
tent with that of the XRB (r ~ 1.4. l Gendreau et all 119951 : 
IChen et all Il997t IVecchi et ail . 119991) . As shown in the pre- 
vious Sections, optically obscured AGN generally have red 
optical color and the average spectral slope obtained for ob- 
sc ured AGN is only s l ightly higher than that of red objects 
of IGeorgakakis etaTl (|2006f ). showing that we are sampling 

similar populations. 

According to IWorslev et all (|2005[ ), whilst the XRB is 
- 85% and ~ 80% resolved in the 0.5 - 2 and 2 - 10 keV 
bands respectively, it is only ~ 60% resolved above ~ 6 keV 
and ~ 50% resolved above ~ 8 keV. The missing popula- 
tion should be made of faint, heavily obscured AGN located 
at redshift of ~ 0.5 — 1.5, and with intrinsic absorption of 



~ 10 23 — 10 24 cm~ 2 . As noted in Section [5j the sources 
classified as SFG do not show any AGN signature in the 
optical and IR. We also find that the fraction of X-ray ab- 
sorbed sources in the SFG class (~ 67%, 16 out of 24) is 
larger than that of X-ray absorbed sources in the type 2 
AGN class (~ 54%, 33 out of 61). Thus sources belonging 
to this class appear to be good candidates to be responsible 
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Table 2. Mean spectral properties from the stacking analysis of optically unobscured, optically obscured AGN and of 
the subclass of SFG objects. N# is fixed to the galactic value. 



for the XRB in the harder X-ray range. We therefore ap- 
plied the stacking analysis to study separately the spectral 
properties of the SFG population in our sample. 

Only 13 SFGs are detected in the pn at off-axis angle 
< 11', so we grouped them in two flux bins. The spectral 
fits of the SFG stacked spectra give T ^ 1.2 — 1.3, with no 
significant differences in the two bins (Table [H and crosses 
in Fig. [T7|) . Therefore the average spectra of the SFGs are 
harder than those of optically obscured AGN (type 2 + 
SFGs) and even harder than the XRB spectrum in the 
same band. If the population responsible for the high en- 
ergy XRB has the same SED properties as the SFG objects 
discussed in this work, they might go unidentified as AGN 
even in the IR, where they look like star forming galaxies. 
A more detail e d disc ussion about this topic is presented in 
iPolletta et ail (j2007fl . 

7. The surface density of optically obscured and 
unobscured AGN 

In Paper I we computed the logN-logS distribution for all 
the sources detected with a probability of false detection 
P < 2 x 10" 5 in the 0.5 - 2 and 2 - 10 keV bands; for 
the 2 — 10 keV band, this probability threshold is slightly 
lower than the 3 a threshold chosen for the present sample, 
therefore all the sources of the 3 a hard sample were in- 
cluded in the logN-logS. We used the differential logN-logS 
for the 2 — 10 keV band reported in Paper I, computed the 
fraction of optically unobscured and obscured AGN in each 
flux bin from the present sample and rescaled the logN-logS 
relationship accordingly. We made the reasonable assump- 
tion that the fraction of optically obscured and unobscured 
AGN should be the same in the area covered by the VVDS 
as well as in the whole XMDS area. 

The differential logN-logS relationships for optically ob- 
scured and unobscured AGN are shown in Fig. [T8l The er- 
rors are the combination of the errors on the original logN- 
logS with those on the fractions, according to the error 
propagation formula. The two logN-logS are quite similar, 
except for the faintest fluxes, where the density of the opti- 
cally unobscured AGN is significantly lower (by a factor of 
~ 4.6) than that of optically obscured AGN. Considering 
the cumulative logN-logS instead of the differential one, we 
can give an estimate of the integral surface density of op- 
tically obscured and unobscured AGN at F > 10~ 14 erg 
cm~ 2 s _1 . We find 138 and 59 sources deg -2 , respectively, 
and the ratio between optically obscured and unobscured 
AGN is ~ 2.3 for the whole flux range covered. The ratio 
would decrease from 2.3 to 1.1 if we assume that the frac- 
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Fig. 18. Differential logN-logS for unobscured (empty dia- 
monds) and obscured AGN (filled diamonds). The units of 
n(F) are number per 10 -15 erg cm -2 s -1 deg -2 . Vertical 
bars show the error on the number of sources, while the 
horizontal bars show the flux bin width. 



tion of unobscured AGN should be corrected by a factor of 
1.6 (see subsection 15. lj) . 

We compared these values with the surface densities 
of broad line and non broad line AGN, estimated by 
iBauer et al.l (|2004h in their study of the Chandra Deep 
Fields. We obtained values from their Fig. 4 and 8. At 
^2-10 > 10 -14 erg cm -2 s _1 , the surface densities differ 
by a factor of ~ 2. This discrepancy is a consequence of the 
fact th at the logN-logS com puted in Paper I is lower than 
that of I Bauer et al.l $2004) . As discussed in Paper I, the 
XMD S log N-logS is slight l y lowe r than those of lBaldi et all 
(|2002fl and lMoretti et alj (|2003f ) for F 2 - W > 2 x 10" 14 erg 
cm" 2 s -1 , but co nsistent within the errors. Instead, the 
IBauer et ail (|2004l ) surface density for fluxes F 2 - 8 ~ 1CT 14 
erg cm" 2 s" 1 is sl ightly higher than that obtained by 
iMoretti et al.l (|2003h . Since all these surveys refer to small 
connected areas in different parts of the sky, it is possible 
that the differences in the derived number counts may be 
due to cosmic variance. 
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Fig. 19. X-ray (2 — 10 keV) absorption corrected lumi- 
nosity vs intrinsic column density. Empty diamonds are 
optically unobscured AGN and filled diamonds are op- 
tically obscured AGN. Encircled points are objects with 
Fx/Fr > 40. The dashed lines mark the region where type 
2 QSOs should be found. 



On the other hand, the iBauer et al. (2004) ratio of non 
broad line and broad line AGN at F ~ 10~ 14 is ~ 1.5 - 2, 
intermediate between our values of 2.3 and 1.1. These au- 
thors describe some caveats about their AGN classification 
criteria, and point out that their number counts of broad 
line AGN must be considered a lower limit and that of non 
broad line AGN an upper limit, so that their non broad 
line/broad line AGN ratio could approach our lower esti- 
mate. 



7.1. Type 2 QSO candidates 

We searched for type 2 QSO candidates in the 3 a hard 
sample. In the X-ray domain the type 2 QSO popula- 
tion is characterized by high intrinsic absorption (N# > 
10 22 cm -2 ) and high X-ray luminosity (L x > 10 44 erg 
s" 1 ). Since locally (in the Seyfert luminosity regime) X- 
ray absorbe d AGN are 4 times more numerous than unab - 
sorbed ones (|Maiolino fe Riekd . 119951 : iRisaliti et "ail 119991) . 
according to the unified AGN model, one would expect 
that the same should be true at high luminosities and 
redshifts, i.e. in the QSO regime. A still undiscovered 
large population of obscured AGN is indee d predicted 
by X - ray background synthesis models (e.g. iGilli et all 
20011: IFranceschini et all l2002t [Gandhi L Fabianl . 120031: 
Uedaet all . 120031 : IWorslev et all I2005D . Before the ad- 
vent of Chandra and XMM - Ne wton, only a fe w type 
2 QSOs were known (see e.g. lAkivama et all . Il998t 



iDella Ceca et all . 120001: IFranceschini et all teOOoTT Peep X 
ray surveys found in deed a large f raction of the objects 
to be obscured (e.g. iBarger et all l2003t iMainieri et all 
2002), but the number of identified type 2 QSOs is still 
very low compared to the model predictions. A signifi- 



cant fraction of high - z, obscured QSOs may have es- 
caped optical spectroscopic identification due to the weak- 
ness of their optical counterparts and misclassification due 
to the lack of AGN signature. On the other hand, medium 
deep X-ray surveys, covering relatively large sky areas at 
a higher flux limit, proved to be effective to select signif- 
icant samples of type 2 QSO candidates among objects 
with hi gh values of the X -ray to optical ratio (Fx /Fr > 
10, see iFiore et all |2003[ ). Recent findings also suggest a 
connection between Extremely Red Objects (EROs, R — 
K > 5 in the Vega system) and type 2 QSOs (see e.g. 



iBrusa et~aL . 20051: Severgnini et all 120051 and references 
therein). iMaiolino et al.l (|2006( ) suggest that by selecting 
extreme values of Fx/Fr(> 40) and extreme values of 
R — K(> 6), the type 2 QSO selection efficiency may ap- 
proach 100%. 

The absorption corrected X-ray luminosity is shown as 
a function of the intrinsic column density in Fig. 1191 There 
is a significant number of objects (34 out of 124) having 
25 L x > 10 44 erg s" 1 and N# > 10 22 cm" 2 in our sample. 12 
of them (35%) are classified as unobscured AGN, while the 
remaining 22 are classified as obscured AGN, based on the 
SED classification. We verified that these objects generally 
have high X-ray to optical ratios, in particular 21 of the 
25 X-ray sources in the 3 a sample having Fx/Fr > 10 
have N/f > 10 22 cm" 2 and L 2 -io > 10 44 erg s _1 . On the 
other hand, while for the optically unobscured AGN the 
fraction of high luminosity, X-ray absorbed sources having 
Fx/F R > 10 is only 25% (3 out of 12), this fraction is 77% 
(17 out of 22) for the obscured ones. 

All the 8 objects with F X /F R > 40 (the threshold 
used by IMaiolino et all l2006f ) satisfy the X-ray definition 
of a type 2 QSO (see encircled objects in Fig. [19]). All 
are optically obscured. These results confirm that type 2 
QSO candidates are found between the high X-ray to op- 
tical ratio population and that the threshold proposed by 
IMaiolino et al.l (|2006l ) is highly efficient in finding type 2 
QSOs but it is far from exhaustive (i.e. many type 2 QSOs 
have F x / Fr < 40). 

In Fig.ED]the color between VVDS R band and SWIRE 
4.5 /jm band is shown as a function of redshift for the type 2 
QSO candidates. Objects fitted by a type 1 AGN template 
have generally blue colors. About 70% of the candidates fit- 
ted by a type 2 AGN or a SFG template have extremely red 
infrared/optical flux ratios, as observed in extremely ob- 
scured AGN and similar to those observed in spectroscopi- 
call y confirmed type 2 QSOs at high redsh i ft (z = 1.5 — 2.5, 
see [Severgnini et all 120061 : iPolletta et all . l2006j ). We have 
K magnitudes from the UKIDSS Early Data Release or 
from the VVDS for 18 objects: 4 of them are EROs, all 
having F X /F R > 40 and all fitted by a type 2 AGN or a 
SFG template. Given the blue optical/IR colors of the op- 
tically unobscured AGN, we exclude them from the type 
2 QSO candidates. Therefore the sample of type 2 QSOs 
reduces to 22 objects. One of them (XMDS 55) has been 
indeed spectroscopically confirmed as a type 2 object (see 
Garcet et al. in prep.). The type 2 QSO candidates repre- 
sent (18 ± 4)% of the sources in the 3 a hard sample and 
have X-ray fluxes in the range 1 — 5 x 10" 14 erg cm" 2 s" 1 . 

As done before for optically obscured and unobscured 
AGN, we rescaled the surface density of X-ray sources at 
F 2 - 10 > 10" 14 erg cm" 2 s" 1 to the type 2 QSO frac- 
tion to estimate the type 2 QSOs surface density. It re- 
sults (35 ± 8) deg" 2 , lower but consistent, within errors, 
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Fig. 20. Optical-infrared color (in the AB system) as a 
function of redshift for the type 2 QSO candidates com- 
pared with the expected colors for various types of galaxy 
and AGN templates. Filled circles are sources fitted by a 
type 1 AGN template, squares are sources fitted by a type 
2 AGN template, asterisks are object fitted by a SFG tem- 
plate. Encircled points are objects having Fx I Fr > 40 
Triang le is SWIRE_J104409.95+585224.8 inlPolletta et al 

2006J), diamond is XBS J0216-0435 in ISevergnini eTal 

2006). 



with (45 ± 15) deg" 2 found bv lCocchia et all (|2007t l in the 
HELLAS2XMM at the same flux level. 

The type 2 QSO population represents about 35% of all 
high luminosity sources (Lx > 10 erg s~ 4 ) in our sam- 
ple. For comparison, according to iPerola et alj (|2004h . in 
the HELLAS2XMM the fraction of X-ray absorbed sources 
(Nfj > 10 22 cm~ 2 ) in the high luminosity (L x > 10 44 erg 
s~ 4 ) AGN population would be between 28% and 40%. 



8. Discussion 

8.1. Optical Obscuration vs X-ray Absorption 

The wide and well sampled multiwavelength coverage from 
the optical through the mid-lR allowed us to use the photo- 
metric approach to redshift determination and classification 
in a very effective way. The comparison with a spectroscopic 
sample gives us confidence in the estimated redshifts and 
in the fact that a photometric type 1 classification is unam- 
biguous, but reveals a bias against the recognition of a num- 
ber of broad line objects as type 1 SEDs. This is due to the 
coexistence of Seyfert 1.8 type SEDs with the presence of 
broad line emission. Since we adopt here the SED classifica- 
tion, these objects will be considered "optically obscured" . 
In optically obscured objects the optical-UV emission from 
the AGN may be either dimmed by intervening dust or be 
weaker than t h at of t he host. In t h is sen s e, "red quasars" 
(I Wilkes et all. l2002t [Gregg et all l2002t iGlikman et al l . 
12004 lUrrutia et all 120051 : I Wilkes et all I2005D would be 
classified here as obscured objects. 

The X-ray data give us an independent and comple- 
mentary information essential to identify AGN where the 
nuclear X-ray emission is heavily absorbed, and AGN fea- 
tures may be completely hidden both in the optical and IR 



bands, as in the case of sources identified with SFGs. In the 
following we discuss the trends of optical obscuration and 
X-ray absorption within our sample separately, in order to 
explore to what extent the two are associated. 

The sample contains 39 optically unobscured AGN and 
83 optically obscured AGN (of which 22 best fitted by a 
SFG template). The two sources with X-ray to optical ra- 
tios and luminosities typical of normal galaxies are excluded 
from this analysis. The X-ray absorbed AGN are 60. X-ray 
absorption occurs in 48 of the 83 obscured AGN (of which 
16 are SFGs). The remaining 12 X-ray absorbed AGN are 
classified as unobscured on the basis of their SEDs. We 
conclude that X-ray absorption is commonly but not ex- 
clusively associated with obscuration since 30% of the un- 
obscured AGN are X-ray absorbed. 

The numbers of optically unobscured, optically ob- 
scured and X-ray absorbed AGN in different flux, redshift 
and luminosity bins are given in Table [31 where we also 
report separately the number of X-ray absorbed sources 
within the optically unobscured and obscured AGN respec- 
tively. 

The fractions of obscured/absorbed AGN in our sample 
are shown in Fig. [^l] as a function of observed flux, redshift 
and absorption corrected hard X-ray luminosity. The left 
panels refer to optically obscured and unobscured AGN, 
while the right panels refer to X-ray absorbed (Wff tr > 10 22 
cm~ 2 ) and unabsorbed AGN, irrespectively of their SED 
classification. 

We used the Bayesian statistics to estimate the "true 
value" o f the fractions and th eir errors (68% credible inter- 
val, see lAndreon et al.l . 120061 , and references therein) and 
to evaluate the reliability of the suggested correlations. We 
tested whether existing data support a model in which the 
fraction 1) is constant or 2) has a linear dependence with 
redshift, flux or luminosity, by computing the Bayesian 
Information Criterion (BIC, Schwartz 1978; an astronomi- 
cal introduction to it can be found in Liddle 2004) for both 
models and then the difference ABIC between the BICs of 
the two models. A ABIC of 6 or more can be used to reject 
the model with the largest value of BIC, wherea s a value 
between 2 and 6 is only suggestive (|Jeffrevsl . 119611) . We also 
compared the trends of optically obscured and X-ray ab- 
sorbed AGN using the same criterion. 

Obscured sources are dominant in the lowest flux bin 
(79lg'Jg%, upper left panel), although a systematic trend 
of optically obscured AGN to decrease with X-ray flux is not 
established (ABIC = 0.7). This result was already appar- 
ent from the logN-logS curves, where the surface density of 
optically obscured AGN largely exceeds that of unobscured 
objects in the lowest flux interval (F < 2 x 10~ 14 erg cm~ 2 
s _1 ), see Fig.[T5J The fraction of X-ray absorbed AGN (up- 
per right panel) is also higher at the faintest fluxes and 
there is a positive indication that it changes systematically 
with flux (ABIC = 4.9). 

The fraction of optically obscured AGN shows instead a 
steep decrease as a function of redshift (middle left panel) , 
from ~ 90% at z < 0.5 to ~ 30% at z > 2, and the 
trend is highly significant (ABIC = 26.3). Similarly, the 
data strongly support a decrease of the fraction of optically 
obscured AGN with luminosity (lower left panel, ABIC = 
23.9). The two trends are not independent, since in flux 
limited surveys a correlation between luminosity and red- 
shift is expected. For X-ray absorbed AGN, data suggest 
constancy with both redshift and luminosity. In summary 
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— 14 < logi^-io < — 13.75 


27 


5 


22 


17 


2 


15 


1 O *7 f -* 1 171 - 1 O f 

— 13.75 < logF2_io < —13.5 


61 


20 


41 


32 


8 


24 


i n C , 1 , ,771 ^- lOOtt 

— 13.5 < log-f^-io < —13.25 


25 


11 


14 


9 


2 


1 


— 13.25 < logF2-io < —12.5 


9 


3 


6 


2 





2 


z < 0.5 


28 


2 


26 


13 





13 


0.5 < z < 1 


42 


11 


31 


14 


1 


13 


1 < z < 1.5 


21 


5 


16 


15 


3 


12 


1.5 < z < 2 


11 


6 


5 


6 


1 


5 








K 
O 


1 9 


7 


o 


logL 2 -io < 42.5 


5 





5 


2 





2 


42.5 < logia-io < 43.5 


27 


1 


26 


15 





15 


43.5 < logL 2 -io < 44.5 


54 


15 


39 


23 


4 


19 


logL 2 -io > 44.5 


36 


23 


13 


20 


8 


12 



Table 3. Total number of sources, number of optically unobscured and optically obscured AGN, number of X-ray 
absorbed AGN (N# > 10 22 cm -2 ), and number of optically unobscured and obscured AGN among X-ray absorbed AGN 
in each X-ray flux, redshift and absorption corrected X-ray luminosity bin (see text and Fig. I21|) . Fluxes are in erg cm -2 
s _1 , luminosities are in erg s _1 . 



there is evidence that the trends of optically obscured and 
X-ray absorbed AGN are different both as a function of 
redshift (ABIC = 17.5) and luminosity (ABIC = 20.1), the 
former showing a decrease with redshift and luminosity, the 
latter being essentially constant. 



8.2. Comparison with the literature 

Several authors compute the fraction of X-ray absorbed 
or optically obscured AGN as a function of all or some of 
the quantities described above, however in the literature a 
detailed comparison between optical obscuration and X-ray 
absorption seems to be lacking. 

The trends of broad line AGN as a function of 
redshift and X-ra y luminosity are e xplored e.g by 
ISteffen et all (l2003f ) . iBarger et al l (|2005h . iTreister &; Urrvl 
(|2005f) . lTozzi et al.l (|2006h . who use data from the Chandra 
Deep Fields, in some cases complemented by shallower 
Chandra observations. All these samples reach flux levels 
significantly deeper than ours thus probing the AGN pop- 
ulation in more depth; given this, we limit ourselves to a 
qualitative comparison. Assuming that unobscured AGN 
correspond to broad line AGN, the trends presented by the 
above authors are in agreement with those shown in the 
middle and lower left panels of Fig. [2TJ 

We examined i n more detail the data of the 
HELLAS2XMM ldf (jFiore et al.l . 120031 ; IPerola et al.l . |2004[ ) 
and of the Serendi pitous Extragalacti c X-ray Source 
Identification (SEXSI. Hckart et all 120061 ). whose flux lim- 
its and areas are comparable to those of the XMDS. For the 
HELLAS2XMM, we computed the fraction of optically ob- 
scured AGN using all the sources for which a spectroscopic 
classification is available (their samples SI, S2 and S4). 
Consistently with our classification scheme, we grouped 
together the objects spectroscopically classified as type 2 
AGN, emission line galaxies (ELGs) and early type galax- 
ies (ETGs), considering them as optically obscured AGN. 
We did the same for the SESXI data, considering as opti- 
cally obscured AGN all sources spectroscopically classified 
as narrow line AGN (NLAGN), ELGs and absorption line 



galaxies (ALGs). Broad line AGN are instead classified as 
optically unobscured. 

We computed the fractions in the HELLAS2XMM and 
SEXSI using the Bayesian statistics, as done for our sam- 
ple. We find that also in the HELLAS2XMM and SEXSI 
cases the fraction of optically obscured AGN decreases (and 
conversely the fraction of optically unobscured AGN in- 
creases) with redshift and luminosity. We notice that in the 
HELLAS2XMM survey the fraction of unobscured AGN is 
larger by a factor of ~ 3 than ours, for redshifts z < 1.5, 
while it is consistent with ours at higher redshifts. The 
agreement with the SEXSI survey is instead better. The 
spectroscopic completeness is however about 90% for the 
HELLAS2XMM, while it ranges from 40% to 70% for the 
SEXSI. The larger fractions of type 1 AGN in spectroscopic 
samples are expected, given the differences between pho- 
tometric and spectroscopic classifications discussed above, 
however the fraction of type 1 AGN in the HELLAS2XMM 
is still larger than ours even when the correction factor com- 
puted in subsection 15. II is taken into account. Nevertheless, 
the trends observed in spectroscopic samples are consistent 
with ours. 

The X-ray properties (fraction of AGN with N H > 10 22 
cm~ 2 ) are explored by a number of authors who use XMM 
Newton or Chandra data of si milar depth (flux limit 



of ~ 10" 



erg cm 



see e.g. 



Piconcelli et al 



IPerola et all 12004 lEckart et all I2006t lAkvlas et al 



and a quantitative comparison with their results is possi- 
ble. Taking into account the different selection criteria and 
corrections that the different authors apply to the data, we 
find agreement within the errors with the results reported 
here. 

Again, we concentrated in particular on the results ob- 
taine d in the HEL L AS2X MM and in the SEXSI su rveys. 
Both IPerola et all (|2004D and lEckart et all (|2006l ) show 
that the fraction of X-ray absorbed AGN increases with 
decreasing X-ray flux, even if the trends are significant 
only when X-ray fluxes as faint as 10~ 15 erg cm -2 s _1 
are considered. In our flux range, we are consistent with 
the HELLAS2XMM and SEXSI values. IPerola et all (2004) 
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and lEckart et al.l (2006) also find that there is no evidence 
of a dependence of the fraction of X-ray absorbed AGN on 
luminosity. Again, these results are consistent with ours, 
with a better quantitative agreement with the SEXSI than 
with the HELLAS2XMM survey (for example, the fraction 
of X-ray absorbed AGN is between ~ 0.4 and ~ 0.6 in our 
case and in the SEXSI, while it is between ~ 0.2 and ~ 0.4 
in the HELL AS2XMM) . 

In conclusion, the analysis of our, the HELLAS2XMM 
and the SEXSI data indicates that the percentages of op- 
tically obscured and X-ray absorbed AGN within the same 
sample show different dependences on redshift and X-ray 
luminosity. Recent models which describe the cosmolog- 
ical evolution of the AG N space density, such as those 
of lUeda et all (|2003l ) and lLa Franca et all (|2005l ). predict 
that the fraction of X-ray absorbed AGN decreases with 
luminosity, and i ncrea se with redshift, in the case of the 
lLa Franca et alJ (|2005l ) model. The combination of a de- 
crease in luminosity and an increase with redshift within a 
single flux limited sample, where the redshift and luminos- 
ity dependences tend to compensate each other, may well 
be the explanation underlying the observed "constancy" 
with redshift and lu minosity of the absorbe d AGN fraction 
in our data. In fact. lLa Franca et alJ (|2005h point out that 
the result of the opposite trends in Lx and z leads to an 
apparent constancy of the X-ray absorbed fraction of AGN. 
Only by combining several samples and thus covering wide 
strips of the Lx — z plane with almost constant redshift 
or luminosity is it possible to disentangle the true depen- 
dences. 

On the other hand, in a recent analysis of the XMM 
- Newton observation of the Chandra Deep Field South, 
iDwellv k. Page] ([20061 ) do not find any dependence of the X- 
ray absorbed AGN fraction on redshift and luminosity and 
suggest that the trends observed by other authors could be 
the result of using deep X-ray data from Chandra, which 
could be biased against high redshift X-ray absorbed AGN. 
Therefore the redshift and luminosity dependence of X-ray 
absorption in the AGN population is still an open issue. 

In any case, the different redshift and luminosity de- 
pendences observed for optically obscured and X-ray ab- 
sorbed AGN imply that in a significant number of objects 
obscuration and absorption are not strictly related, more- 
over the relation depends on redshift/luminosity in a sys- 
tematic way. 

8.3. Objects with different absorption properties 

There is a number of examples in the literature of ob- 
jects that have opposite X-ray and op tical properties , 
such as X-ray absorbed type 1 AGN (e.g . Comastri et all 
120011: iBrusa et all . 120031 : lAkivama et all 120031) or X-ray 



ma et an izuudi) or A-ray 
Panessa fc Bassaui. l2002t 



unabsorbed type 2 AGN (e.g. ^ . . _ _ _ 

ICaccianiga et all 12004 IWolter et all l2005f) . however it is 
not clear so far how common these exceptions are and how 
they c an be reconci l ed wit h the unified model (jAntonuccil . 
119931) . iPerola et alJ (|2004l ) find tha t about 10% of br oad 
line AGN are X-ray absorbed, while iTozzi et alJ (|2006| ) es- 
timate that the correspondence of unabsorbed (absorbed) 
X-ray sources to optical type 1 (type 2) AGN is accurate for 
at least 80% of the sources. We address this question in the 
following subsections, where instead of type 1 and type 2 
AGN as derived from the optical spectra, we will consider 



optically obscured and unobscured AGN as derived from 
the SED classification. 



8.3.1. X-ray absorption without optical obscuration 

We find 12 X-ray absorbed, optically unobscured AGN, 
31% of all unobscured AGN. In 7 cases, the fit requires 
additional extinction, Ay=0.40-0.55, but lower than that 
would be derived from the column density in the X-ray 
spectral fits assuming the standard dust-to-gas ratio. The 
discrepancy between optical and X-ray properties can be 
explai ned e.g. by a dust-to-g as ratio lower than the Galactic 
value (|Maiolino et all 1200 ll ) or by a different path of the 
line of sights to the X-ray and the optical sources (e.g. see 
dual absorber model in Risaliti et al. 2000 or torus clumpy 
model, Hoenig et al. 2006). Another possibility is that in 
some objects the absorbing gas is ionized rather than neu- 
tral: in that case, dust would likely be absent and the intrin- 
sic continuum plus broad emission lines would be observed 
in the optical spectrum of the AGN. Examples of broad 
line AGN whose X-ray spectra show absorptio n well fitted 
by an ionized absorber model are reported in iPage et all 
(200j). 



8.3.2. Optical obscuration without X-ray absorption 

35 out of 83 optically obscured sources do not show high 
X-ray absorption. We note however that 22 of them (63%) 
are indeed fitted by the Seyfert 1.8 template, which in a 
number of cases corresponds to objects with broad emission 
lines (subsection I5.1[) . It is therefore likely that an interme- 
diate class between unobscured and truly obscured AGN 
exists, made of objects which are dominated by star-light 
emission in the near-IR, and X-ray unabsorbed. In a future 
work we will extend our analysis to the whole XMDS area; 
the larger statistics will allow us to refine the photomet- 
ric classification using a wider set of templates and better 
explore the relation between X-ray absorption and the opti- 
cal to mid-IR SED. However, we do not expect all obscured 
AGN that are unabsorbed in the X-rays to be misclassified. 

Another plausible explanation for (apparent) obscura- 
tion without X-ray absorption is a larger relative luminosity 
of the host galaxy compared to the AGN optical light. In 
this scenario, the AGN optical light might simply be fainter 
than that of the host galaxy and not extincted. 

Our analysis shows that about 40% of our sources have 
opposite optical and X-ray properties (12 X-ray absorbed, 
optically unobscured AGN and 35 X-ray unabsorbed, op- 
tically obscured AGN). The uncertainties related to this 
fraction are linked, on the one side, to the large errors in- 
volved in the computation of intrinsic column densities for 
high redshift AGN (see subsection 16. 2[) . and on the other 
side, to our classification of sources based on SED fitting 
templates. However, these effects cannot fully account for 
the large number of objects with discrepant optical and X- 
ray properties and the very different trends that we observe 
in Fig. [2D We remind that similar discrepancies are also ap- 
parent in the literature, where spectroscopic classifications 
are used (subsection 18. 2\ . These results suggest that the 
basic formulation of the unified model, in which the view- 
ing angle is the sole factor in determin ing the AGN type, 
might be too simplistic. As an example, Elitzur (2006) pro- 
pose that the difference between type 1 and type 2 AGN is 
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instead an issue of probability for direct view of the AGN 
through a clumpy, soft - edge torus. Moreover, he suggests 
that the "X-ray torus" does not coincide with the "dusty 
torus" and that the bulk of the X-ray absorption likely 
comes in most cases from clouds located in the inner, dust 
free portion of the X-ray torus. The trends with luminosity 
described above agree with this kind of picture, where at 
high luminosities dust can be evaporated/expelled, while 
absorption by gas can be associated with strong outflowing 
winds outside the dusty torus. 

Larger samples obtained by combining sub-samples 
from surveys of various depths and areas, joined with opti- 
cal spectroscopic data, are necessary to minimize selection 
effects, provide a deeper understanding of the AGN prop- 
erties and test this scenario. 

9. Summary and conclusions 

We have selected 136 X-ray sources detected at > 3<r in 
the 2 — 10 keV band in the XMDS area also covered by the 
VVDS. More than 90% of the sources have been identified 
with an optical and/or infrared counterpart (Section[3]) and 
98% of them have X-ray to optical and X-ray to infrared 
ratios typical of AGN (Section [3]). 

We used the optical and infrared data to construct SEDs 
and compute photometric redshifts. The comparison with 
the spectroscopic redshifts available shows that in 90% of 
cases there is agreement between photometric and spec- 
troscopic estimates (Section[5|). All sources fitted by a SFG 
template which have X-ray to optical ratios typical of AGN 
(22 out of 24) have high hard X-ray luminosities (> 10 42 
erg s _1 ), suggesting that all are indeed AGN with the 
host galaxy emission dominating in the optical - IR bands. 
Objects fitted by a type 1 AGN template have generally 
blue optical/IR colors and in most cases do not show X-ray 
absorption, while those fitted by a type 2 AGN or a SFG 
template have red optical/IR colors and most of them are 
X-ray absorbed (Section [5] and [6]). 

Comparison between photometric and spectroscopic 
classification, when available, shows that the type 1 AGN 
photometric classification is unambiguous, but we under- 
estimate the fraction of broad line AGN, since 10 out of 
26 are classified as type 2 AGN due to the dominance of 
star light emission in the near-IR. AGN fitted by a type 1 
template are referred as optically unobscured, while those 
fitted by a type 2 AGN or a SFG template are referred as 
optically obscured (subsection l5.1[) . 

We extracted the X-ray spectra of the 55 X-ray sources 
having at least 50 net counts in the 2 — 10 keV band. For 
sources with a smaller number of detected counts, we used 
hardness ratios to compute the column densities. We find 
that, when the redshift dependence is taken into account, 
60% of the optically obscured AGN are X-ray absorbed, but 
also 30% of the optically unobscured AGN have N^ tr > 
10 22 cm -2 , showing that optical and X-ray classifications 
are not strictly related (Section O . 

We constructed stacked X-ray spectra to measure av- 
erage spectral properties of our sample and to find differ- 
ences between optically obscured and unobscured AGN as a 
function of X-ray flux. We find that stacked spectra of opti- 
cally unobscured AGN have a photon index consistent with 
r = 2, similar to average values found for X-ray unabsorbed 
AGN. On the other hand, the slope of stacked spectra of 
optically obscured AGN is consistent with T ~ 1.6. The 



stacked spectrum of the objects fitted by a SFG template 
is even harder, with T ~ 1.2 — 1.3. (subsection 16. 3[) . 

Comparing the fractions of optically obscured and X- 
ray absorbed AGN, we find that while the fraction of op- 
tically obscured AGN steeply decreases with redshift and 
luminosity, that of X-ray absorbed AGN is nearly constant 
at ~ 50% as a function of redshift and X-ray luminosity 
(Section [S]). The constancy of the population of X-ray ab- 
sorbed AGN with redshift an d luminosity obse r ved in our 
sample can be explained by the lLa Franca et all (|2005f ) pre- 
dictions that the fraction of X-ray absorbed AGN should 
decrease with luminosity and increase with redshift, since 
these two dependences tend to compensate each other in a 
single, flux limited sample. The different trends of optically 
obscured and X-ray absorbed AGN are confirmed also by 
the analysis of spectroscopic samples from the literature, 
showing that this result is not biased by the uncertainties 
in the photometric classification. 

In 39% sources (47 out of 122) an inconsistency between 
X-ray absorption and optical obscuration is observed (12 
X-ray absorbed, optically unobscured AGN, and 35 X-ray 
unabsorbed, optically obscured AGN). About 63% of the 
optically obscured X-ray absorbed AGN can be indeed mis- 
classified broad line AGN. On the other hand, the signifi- 
cant fraction of optically unobscured, X-ray absorbed AGN 
found in our sample suggests that the basic formulation of 
the AGN unification model can be too simplistic. 

We also computed the differential logN-logS relationship 
for the obscured and unobscured AGN, finding that the op- 
tically obscured AGN begin to dominate for F < 2 x 10~ 14 
erg cm" 2 s , where the ratio between obscured and unob- 
scured AGN is ~ 4.6 (Section [7]). In the whole flux range 
considered, the surface density of the optically obscured 
AGN is higher than that of the optically unobscured ones 
by a factor of ~ 2.4. However, if a correction is applied 
to account for misclassified type 1 AGN, the ratio between 
optically obscured and unobscured AGN is ~ 1.1. 

We find 22 sources that could be classified as type 2 
QSO candidates (N H > 10 22 cm- 2 , L x > 10 44 erg s" 1 ). 
They are fitted by a type 2 AGN or a SFG template and 
on average their infrared/optical and X-ray/optical flux ra- 
tios are typical of extremely obscured AGN. 4 of the 18 
having a measured K magnitude arc EROs (R — K > 5 in 
the Vega system). We estimate a surface density of type 2 
QSOs at F 2 _io > 10~ 14 erg cm- 2 s" 1 of (35 ± 8) deg" 2 
(subsection I7.ip . 

In this work the full power of multiwavelength obser- 
vations is exploited to understand the global properties of 
AGN. We plan to extend our analysis to the whole XMDS 
to improve the statistical significance of our results. 
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Fig. 21. Left panels: fraction of optically obscured (filled diamonds) and unobscured (empty diamonds) AGN as a 
function of X-ray flux (upper panel), redshift (middle panel) and luminosity (lower panel). Right panels: fraction of 
X-ray absorbed (N# > 10 22 cm~ 2 , hlled triangles) and unabsorbed (empty triangles) AGN as a function of X-ray flux 
(upper panel), redshift (middle panel) and luminosity (lower panel). 
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Appendix A: Properties of the 3 a hard sample 

We report in Table[Othe X -ray, optical and infrared prop- 
erties of the 3 a hard sample. Column 1 list the XMDS 
identifier, column 2 the XMDS name, columns 3 and 4 the 
X-ray coordinates (astrometrically corrected) , column 5 the 
2 — 10 keV flux, column 6 the hardness ratio between the 
2 — 10 and 0.5 — 2 keV bands, column 7 the distance be- 
tween the X-ray and the optical VVDS position, column 8 
the VVDS B magnitude, column 9 the VVDS I magnitude, 
column 10 the color classification (B for blue, i.e. B — I<1 
and R for red, i.e. B - I > 1), column 11 the CFHTLS i' 
magnitude, column 12 the SWIRE 3.6 /im flux, column 13 
the SWIRE 24 /im flux, column 14 the spectroscopic red- 
shift, column 15 the photometric redshift, column 16 the 
photometric classification based on SEDs (see text) and 
column 17 the absorption corrected 2 — 10 keV luminosity. 

A means no data available (i.e. source outside the 
field of view or not measured flux because of instrumental 
problems). For sources undetected at 24 //m, a 5 a upper 
limit (241 /xJy) was used. For sources with ambiguous op- 
tical - IR identification, two or more rows are associated 
with the X-ray source, one for each candidate counterpart, 
with the one having the best probability (see Section [3]) 
listed as first. In these cases, X-ray data are reported only 
in the first row. All magnitudes are in the AB system. The 
systematic error on SWIRE fluxes is 5%. 



Table A.l. The 3 a hard sample: XMDS identifier and name, X-ray coordinates, X-ray flux (2 — 10 kcV), hardness ratio between the 2 — 10 and 0.5 — 2 keV 
bands, distance between the X-ray and optical VVDS positions, VVDS B and I magnitudes, color classification, CFHTLS i' magnitude, SWIRE 3.6 and 24 /im 
fluxes, spectroscopic and photometric redshifts, SED classification (see text) and absorption corrected X-ray luminosity. All magnitudes are in the AB system. 
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J022531.4 


-044210 


02 


25:31.431 


-4 


42:10 


616 


2.71 


1.00 


2 


7 


23.86 


22.06 
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22 


28 


34 
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-045740 


02 
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-4 


57:38 
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4.72 
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9 


22.04 


19.46 


R 


19 


64 


340 
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44.02 
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XMDS 


,1022332.0 


-045740 


02 


23:32.191 


-4 


57:38 


712 


2.45 


-0.56 





6 


20.32 


19.62 


B 


19 


66 


262 


1250 
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AGN2 


44.07 
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XMDS 


J022326.3 


-045708 


02 


23:26.549 


-4 


57:05 


982 


2.91 


-0.31 


1 


5 


20.98 


20.48 


B 


20 


32 


132 


579 


0.826 





839 


AGN2 


43.98 




787 


XMDS 


J022317.9 


-045527 


02 


23:18.081 


-4 


55:25 


250 


1.42 


0.18 





5 


24.22 


22.65 


R 


23 


23 


14 


< 241 







981 


AGN2 


43.90 




788 


XMDS 


J022353.7 


-045510 


02 


23:53.880 


-4 


55:08 


174 


4.86 


-0.52 





4 


22.29 


21.43 


B 


21 


52 


51 


< 241 







958 


AGN1 


44.36 




789 


XMDS 


J022329.1 


-045452 


02 


23:29.341 


-4 


54:50 


778 


2.40 
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1 


2 


20.44 


19.92 


B 


20 


01 


140 


< 241 
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AGN2 


43.56 
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XMDS 
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-045120 


02 
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-4 


51:18 
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2.03 


-0.58 
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22.70 


20.95 


R 


21 


07 
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AGN2 


43.88 
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XMDS 


J022344.4 


-045120 


02 


23:44.634 


-4 


51:18 
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1.74 


-0.42 


1 





22.36 


21.50 


B 


21 


05 


29 
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AGN2 


43.92 
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XMDS 


,1022333.0 


-044924 


02 


23:33.185 


-4 


49:22 


513 


1.78 


-0.49 
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20.26 


B 


20 


22 


37 
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2 
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AGN1 


44.87 


ij 
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02 
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-4 
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1.69 
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B 
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15 
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2 
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44.91 
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J022319.4 


-044732 


02 


23:19.582 


-4 


47:30 
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3.87 


-0.38 





7 


21.18 


18.48 


R 


18 


68 


531 


3756 


0.293 
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AGN2 
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c 

H 
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XMDS 


,1022330.6 


-044633 


02 


23:30.803 


-4 


46:31 
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3.35 


-0.50 





7 


24.57 


23.37 


R 






18 
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1 
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AGN2 


44.84 


S- 


828 


XMDS 


J022318.8 


-044616 


02 


23:19.040 


-4 


46:13 


919 


1.73 


-0.51 





1 


21.23 


20.76 


B 


21 


15 


45 


356 







673 


AGN1 


43.53 




840 


XMDS 


J022330.9 


-044235 


02 


23:31.149 


-4 


42:33 


036 


2.07 


-0.21 
















9 


< 241 




2 


297 


AGN2 


44.94 




842 


XMDS 


J022402.4 


-044140 


02 


24:02.650 


-4 


41:38 


339 


3.63 


0.16 


3 


8 


14.87 


13.54 


R 






13528 


40594 


0.043 





010 


SFG 


41.24 




844 


XMDS 


J022343.2 


-044105 


02 


23:43.453 


-4 


41:03 


352 


1.95 


0.75 


2 


8 


26.04 


24.47 


R 


24 


18 


26 


363 




1 


493 


AGN2 


44.57 


> 
























4 


5 


23.70 


23.24 


B 


23 


25 








1 


970 


AGN2 


44.87 


c 


846 


XMDS 


J022317.2 


-044035 


02 


23:17.415 


-4 


40:33 


523 


3.72 


-0.38 


1 


1 


18.81 


18.23 


B 


18 


74 


580 


4587 


0.842 





765 


AGN1 


44.11 




1197 


XMDS 


,7022720.2 


-045738 


02 


27:20.271 


-4 


57:38 


992 


1.90 


0.32 


1 


5 


23.05 


21.70 


R 


22 


48 


39 


306 




1 


116 


AGN2 


44.18 


£ 


1199 


XMDS 


J022651.6 


-045714 


02 


26:51.736 


-4 


57:14 


797 


6.39 


-0.44 





8 


21.93 


19.78 


R 


20 


11 


104 


830 


0.331 





290 


AGN2 


43.37 


r+ 


1201 


XMDS 


J022723.4 


-045608 


02 


27:23.526 


-4 


56:08 


628 


1.65 


-0.43 


1 


5 


22.18 


21.34 


B 


21 


48 


54 


313 







168 


AGN2 


42.12 


O 


1219 


XMDS 


,1022701.6 


-045158 


02 


27:01.690 


-4 


51:58 


996 


1.15 


0.07 
















23 


< 241 




2 


112 


SFG 


44.62 


x 


1226 


XMDS 


J022711.7 


-045038 


02 


27:11.761 


-4 


50:38 


630 


8.55 


-0.51 


1 


1 


21.17 


20.42 


B 


21 


10 


97 


734 







946 


AGN1 


44.59 




1227 


AiVlDb 


,1022736.8 


-045033 


02 


o *7 . o c one 

27:36.895 


— 4 


50:34 


200 


2.41 


0.00 





9 


21.86 


19.44 


R 


19 


60 


109 


864 







445 


AGN2 


43.26 


c 


1231 


XMDS 


J022731.9 


-044957 


02 


27:31.946 


-4 


49:57 


627 


0.98 


-0.54 





8 


22.63 


20.95 


R 


21 


22 


73 


416 







741 


AGN2 


43.39 


cc 


1236 


XMDS 


J022729.0 


-044857 


02 


27:29.058 


-4 


48:57 


908 


1.18 


0.09 


2 


5 


27.35 


24.03 


R 


24 


71 


30 


621 




1 


513 


AGN2 


44.28 




1246 


XMDS 


J022712.8 


-044636 


02 


27:12.916 


-4 


46:37 


102 


2.34 


-0.59 


1 


5 


18.03 


17.69 


B 


17 


75 


451 


2830 




1 


446 


AGN1 


44.48 




1247 


XMDS 


J022633.1 


-044637 


02 


26:33.213 


-4 


46:38 
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3.17 


-0.32 


2 





25.01 


21.29 


R 


21 


06 


72 
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1 
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AGN2 


44.42 




1248 


XMDS 


J022725.4 


-044619 


02 


27:25.522 


-4 


46:19 


686 


4.62 


-0.54 





1 


17.02 


15.29 


R 


15 


36 


055 
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SFG 


41.09 
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XMDS 


J022716.0 


-044539 


02 


27:16.045 


-4 


45:39 
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4.62 


-0.60 


1 


6 


20.42 


20.02 


B 


20 


12 
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43.82 
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J022751.3 


-044251 


02 


27:51.403 


-4 


42:51 


783 


3.19 


-0.36 


1 


7 


20.19 


19.82 


B 


19 


77 


63 
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1 
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AGN1 


44.79 


& 
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XMDS 


J022712.6 


-044221 


02 


27:12.698 


-4 


42:21 


727 


6.28 


0.10 


1 


8 


19.79 


17.90 


R 


18 


03 


435 


4214 
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232 


SFG 


42.93 


=_ 
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Appendix B: Results of spectral fits 

We report here the spectral fit results from spectra of in- 
dividual sources (see section l(TTj) . In Table |B~T1 we list the 
column density, photon index, reduced x 2 and number of 
degrees of freedom obtained from analysis of all sources for 
which we could leave free both T and Njj, together with 
the hardness ratio values. In Table HOI we report spectral 
fit results obtained by fixing the photon index to T = 2.0 
(columns 2 and 3) and T — 1.7 (columns 4 and 5), or insert- 
ing photometric (or spectroscopic, when available) redshift 
in the model (columns from 6 to 9). In column 10 we report 
the EPIC cameras used to extract the spectrum. 

For three sources (XMDS 161, 282 and 1199, see 
Tabic IB.2P the spectral model with r = 2.0 gives a poor 
fit (xl > 2 ). XMDS 1199 shows a moderate X -ray absorp- 
tion (N# ~ 10 21 cm -2 ) and T frozen to 1.7 gives a better fit 
{xl = 1-35). The other two sources have instead very steep 
spectra, in fact the spectral fit obtained with free photon 
index (see Table IB.1[) gives T > 2 for both of them, but, 
while for XMDS 282 the fit with T free is good (xl = 1.07), 
we were not able to find an acceptable fit with a simple 
power law model for XMDS 161, whose spectrum exhibits 
a significant soft excess. 

On the other hand, when T = 1.7, there are 10 sources 
having xl > 2: all of them are bright, soft sources, for which 
the spectral fit with free photon index gives T > 2. 

For sources XMDS 124 and 779 no stable solutions were 
found fixing the photon index to 2.0. 

When the redshift was introduced in the spectral model, 
the photon index was left free when the % 2 statistics could 
be used, otherwise it was fixed to 2.0. In the two cases 
cited above (XMDS 124 and 779) and for XMDS 739 we 
however had to fix the photon index to T = 1.7, because 
for T = 2.0 no stable solution was found. For XMDS 453 
no stable solution was found with either value of T. 
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XMDS ID logNg (cm' 2 ) F x l {d.o.i.) HR cb HR^T 



112 


20.44+nni 
w — 0.03 


1 


65+ ^ 

^^—0.16 


0.90 


(19) 


-0.60 


-0.65 


120 


21.08+°,',? 


1 


—0.19 


1.02 


(32) 


-0.61 


-0.48 


140 


20.41+°'??, 

— U.OU 


1 


60+°?° 

— 0. Id 


0.84 


(10) 


-0.72 


-0.33 


142 


20.41+°™ 

— U.OU 


2 


90 + 0. 42 
— 0.24 


1.15 


(10) 


-0.74 


-0.37 


144 


20.98+°,' » 

^ ^^ — 0.30 


2 


63+?,?° 

^^—0.26 


1.00 


(12) 


-0.60 


-0.51 


161 


w — 0.00 


2 


yy+0.O8 

' ' —0.05 


1.40 


(74) 


-0.81 


-0.60 


179 


20.41+°'*° 

—0.00 


2 


34+035 

—0.25 


1.30 


(11) 


-0.65 


-0.52 


281 


20.97+°'?? 

w ^ ' —0.35 


2 


oy+0.25 
" —0.47 


0.69 


(12) 


-0.72 


-0.33 


282 


20.41+SSS 

— U.OU 


2 


45+°?^ 

— 0. Id 


1.07 


(19) 


-0.72 


-0.62 


416 


20.41+°« 
— u.uu 


2 


n q+0.36 


1.24 


(13) 


-0.60 


-0.52 


440 


20.68+2't? 


2 


QO+ - 50 

u — 0.29 


0.34 


(10) 


-0.58 


-0.52 


561 


20.42+nnn 

w — 0.00 


2 


46+?'" 

—0.21 


0.74 


(13) 


-0.76 


-0.55 


577 


20.41+°'?,?, 

— U.UU 


2 


oc+0.19 
—0.13 


0.60 


(25) 


-0.74 


-0.52 


738 


20.41 + nS 
w — 0.00 


2 


oe+0.28 
ou — 0.26 


0.43 


(8) 


-0.69 


-0.39 


742 


20.41+°'?,? 

— 0.00 


2 


1 o+0. 38 
—0.32 


1.26 


(7) 


-0.66 


-0.50 


788 


20.4ll°;« 


1 


81 + ' 27 
—0 20 


0.62 


(11) 


-0.56 


-0.80 


789 


90 49+0-54 


2 


09+0:49 
Oz -0.22 


0.58 


(11) 


-0.74 


-0.47 


820 


20 41 +°' 61 


1 


4 o+0.51 
48 -0.22 


1.21 


(8) 


-0.52 


-0.41 


1199 


91 OO+0.23 
0.30 


2 


1 O+0.32 
io -0.32 


1.06 


(9) 


-0.58 


-0.32 


1226 


20.54+°;?° 


1 


OO+0.21 
8 ^-0.15 


0.80 


(22) 


-0.62 


-0.49 


1248 


20 41 +°' 47 


1 


no+0.23 
yo -0.15 


1.01 


(20) 


-0.67 


-0.33 


1252 


90 41 +°' 21 


2 


9n +0.20 
^ u -0.19 


1.24 


(15) 


-0.66 


-0.69 



Table B.l. Spectral parameters obtained using simple power law model with both T and % free parameters for sources 
for which \ 2 statistics can be used. The quoted errors correspond to the 90 per cent confidence level for one interesting 
parameter. The formal 0.00 errors on logN# are the result of having fixed a minimum column density (the galactic value). 
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Results if 2 
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r = 


2.0 


r = 


1.7 
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logN w (cm -2 ; 


i X 2 . (d.o.f.) 


logNn (cm -2 ; 


i x2 (d.o.f.) 


logN H (cm" 2 ) 


r 


X 2 (d.o.f.) 


z 






4 


21.64lS-.i2 


- 


21.52l_? 


- 


2i.8iiS:iS 


2.0 


- 


0.265" 


pn 




18 


21.26± ;i| 


- 


21.07«; 2 » 


- 


00 ■■ 9 +0.19 
— 0.24 


2.0 


- 


1.612 


pn 




36 


21.48l°;„ 


- 


91 *}1 +0-23 
4±.o±_ 32 


- 


00 1 Q+0.25 
^^.±o_ 29 


2.0 


- 


1.043 


Ml 


M2 


40 


ryr\ 

ZZ.ZZ_q.23 


- 


22 14+ ' 24 

ZZ. l^_o .25 


- 


23.42l°; 23 5 


2.0 


- 


1.955 


pn 




55 


21.86lS;S 


- 


91 7 o+0.29 
^1. I O_o.38 


- 




2.0 


- 


1.029° 


pn 




60 


20.42lg;« 


- 


20.41±o ; gi 


- 


20 41 + ' 53 

ZU.fll_ oo 


2.0 


- 


0.478" 


pn 




71 


21.94±g;_ 


- 


21.8ll° 22 


- 


22 91+° 19 


2.0 


- 


1.401 


Ml 


M2 


91 


on 41 + - 23 


- 


20 41+° 16 
zu.^i_ .oo 


- 










pn 




94 


21.77±g;„ 


- 


21.67l 1 ;i 7 


- 


22.03±g;i? 


2.0 


- 


0.349 


pn 




111 


21.50±g;_ 


- 


91 oc+0.25 
Z1,OJ -0.33 


- 










pn 




112 


21.00±g;„ 


1.01 (20) 


20.58+^ 7 


0.86 (20) 


20.61±g;|g 


l-68l°o^ 


0.89 (19) 


0.878" 


Ml 


M2 


120 


21.09±g;iS 


0.99 (33) 


20.75l°i 8 5 


1.16 (33) 


20.73i°;g 








Ml 


124 






20.90+°;^ 


- 


1.7 


- 


0.371 


pn 


M 

a> 


133 


21.50±g;l? 


- 


9 -| O 7 + 0.18 

-0.23 


- 


21.98±g"if 


2.0 


- 


0.723 


pn 


139 


on 49+0 31 


- 


90 41 +° 15 


- 


20 41 + ' 47 

ZU.tl_ go 


2.0 


- 


0.729 


ml 


Q3 

M? 


140 


20.93+034 


1.45 (11) 


20.48+oo7 


0.85 (11) 


on 41 +0 36 
ZU.tl_ oo 


1 60+ - 21 


0.84 (10) 


0.22" 


pn 


> 


142 


20 41+° :2 ° 


1.38 (11) 


20 42+ - 00 


2.75 (11) 


90 41 +°' 66 
ZU.11_ oo 


9 9O+0.37 
2.28_ 24 


1.15 (10) 


0.573 


pn 




144 


on 49+0-27 


1.25 (13) 


2o.4it° :S2 


2.26 (13) 


21-36±g;„ 


9 ci +0.50 
^• u± -0.26 


0.99 (12) 


0.864 


pn 


5' 


149 


20.93+047 




20.5ll°f 




21.081°;?* 


2.0 




0.493 


pn 


ty 


161 


20 42+ o:o ° 


13.24 (75) 


20 4i+° oo 

ZU.'±l_ oo 


22.11 (75) 


20 41+ 00 

ZU.tl_ go 


9 77+0.07 
1 1 -0.05 


1.40 (74) 


0.053" 


pn 


CD 


179 


on 41 +°- 14 

ZU.11_ oo 


1.63 (12) 


20 42+ 00 


2.96 (12) 


20 41 +°' 45 

ZU.tl_ go 


2 04+0.32 

Z -°^-0.25 


1.30 (11) 


0.495" 


pn 


X 


a : spectroscopic redshift 






















Table B.2. Spectral parameters obtained using a simple power law model with z = and T = 2.0 in columns 2 and 
3, r = 1.7 in columns 4 and 5 and using photometric (or spectroscopic when available) redshift in colums 6, 7 and 8. 
In this case, the photon index is a free parameter when \ 2 statistics was used. The quoted errors correspond to the 90 
per cent confidence level for one interesting parameter. The formal 0.00 errors on logN# are the result of having fixed a 
minimum column density (the galactic value). 
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logNff (cm" 2 ; 


) X 2 (d.o.f.) 


logNff (cm -2 ) 


i xl (d.o.f.) 
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r 


xl (d.o.f.) 
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272 


20.41_°; d 


- 


2041 +o.09 
zu.i±_ .o 


- 


20 4i+ u - 1 ' 4 
zu.^i_ oo 


2.0 


- 


0.797 


pn 




281 


20.571?;?? 


0.75 (13) 


90 41 + 18 
ZU.41_ 000 


1.34 (13) 


91 QC-+0.36 
Zl.00_o.66 


9 09+O.I9 
Z.OZ_ 22 


0.67 (12) 


0.93 


pn 




282 


20 41+ - 00 


2.18 (20) 


20 42+ 00 
zu -^ z -o.oo 


4.71 (20) 


20 41 +°' 32 

zu - 41 -0 .00 


9 4 r+0.26 
z -*°-0.16 


1.07 (19) 


0.189 


pn 




406 


oo 9Q+ - 21 

zz.zt/_Q 24 


- 


99 90+° 
zz. zu_o 25 


- 


22.88l°- 2 4 


2.0 


- 


0.713 


pn 




416 


20.421°;^ 


1.16 (14) 


20 41 +° 19 

ZU.41_o.00 


1.55 (14) 


2o.4ii;-.s; 


Z - UO _0.23 


1.24 (13) 


1.668 


pn 




440 


20.651^3 


0.31 (11) 


90 49+ ' 25 
zu -^ z -o.oo 


0.62 (11) 


21.501°;^ 


2.00l°j? 


0.34 910) 


2.574 


pn 




453 


99 9/1+0-39 
zz - zy: -0.38 


- 


22 10+ ' 41 

— 0.66 


- 










M2 


pn 


498 


90 49+ - 44 


- 


90 41 + - 23 
ZU. 41_ oo 


- 


20.4ll ; 2 ° 


2.0 


- 


2.031" 


pn 




521 


91 4Q+ ' 22 


- 


91 Qtr+0.26 
Zl.OO_ 4 6 


- 


99 oq+0.17 
zz.z» _ 40 


2.0 


- 


1.325" 


Ml 


M2 


561 


20 41+° 12 
zu.ii_ o 


1.61 (14) 


20 41 + ' 00 
ZU. 11_ .00 


3.26 (14) 


20.42l ;g 2 


2.48l°;g 


0.74 (13) 


1.678 


pn 




564 


20.42+°-^ 


- 


20.4ll '; 4 


- 


20.4llH 3 o 


2.0 


- 


0.224 


Ml 




577 


90 4i+° : ' )9 


1.31 (26) 


20.411°°° 


3.08 (26) 


20.411°;^ 


9 qc+0.18 
z - od -0.13 


0.60 (25) 


0.252 a 


Ml 


MJ 

(t> 
M 

o 

c+ 


710 


20.83.__2 


- 


90 4.9+ - 26 
zu -^ z -o.oo 


- 


21.741?; 24 


2.0 


- 


2.320 


pn 


718 


21.29l°; 27 , 


- 


21 09+°' 34 


- 


21.811°;^ 


2.0 


- 


0.896 


pn 


738 


20.42lg;« 


0.97 (9) 


20 41 +°' 00 

zu ' 41 -0 .00 


2.40 (9) 


20.41_°;°_ 


9 Qfi+0.44 
z.ou_ 26 


0.43 (8) 


1.010" 


pn 


739 


22.31_°;__ 


- 


21 Q4+ - 64 

ZJ -- 34: -0.52 


- 


22.051°,;.. 


1.7 


- 


0.140" 


pn 


> 


742 


90 41 + - 28 


1.16 (8) 


90 41 + - 21 
zu.ii_ 00 


1.75 (8) 


20 41 + 1 06 

ZU.11_o oo 


9 i o+0.38 


1.26 (7) 


1.924" 


pn 


a 


746 


in 1 a 

22.50l°;i* 




99 41 +° 15 

— 0. 16 




99 71 +0-15 
zz - ' i -0.15 


2.0 




0.219 


Ml 




755 


20.4ll°; 2 « 




90 41 + 15 




20 41+ ' 41 

ZU.11_ oo 


2.0 




0.730 


pn 


ty 


779 






91 t^S+039 
Zl.OO_ 3 




22.10l°i° 


1.7 




0.637 


pn 


CD 


782 


20.421°;^ 




20 42+ ' 32 
zu -^ z -o.oo 




20.411°;™ 


2.0 




0.826" 


pn 


X 



a : spectroscopic redshift t_ 



Table B.2. continued 

O 

s 



XMDS ID 




Results if z = 








Results if 2 


Zph 






EPIC" 




r = 


2.0 


r = 


1.7 












camc^ 




logN_- (cm~ 2 ~ 


I xl (d.o.f.) 


logN ff (cm" 2 ) 


i xl (d.o.f.) 


logNtf (cm" 2 ) 


r 


Xl (d.o.f.) 


Zph 




788 


20.61_S;£ 


0.71 (12) 


90 41 + u - 2 ' 
ZU.11_o. o 


0.63 


(12) 


90 41 + u - y 3 
zu.ii_ oo 


1 01 +0.24 
1 -° i -0.20 


0.62 


(11) 


0.958 


pn g- 


789 


20.4ll°; 6 


1.00 (12) 


20 42+°-°° 
zu.iz_ oo 


2.45 


(12) 


20 41 +°- 76 
zu.ii_ go 


9 oo+0. 41 
z - 00 -0.23 


0.58 


(11) 


0.604" 


pn ^" 


820 


21.08l°; 22 


1.41 (9) 


20.79l°'3? 


1.14 


(9) 


20 41+ ' 74 
zu.ii_ oo 


l.OU_ 23 


1.21 


(8) 


0.293" 


pn Qd 


825 


21 41+° 12 

Z1.11_ 14 




91 99+O.I6 
Zl.zz_ 21 






22.36+°-^ 


2.0 






1.735 


Ml M2 


828 


20 42+ ' 41 

zu -^ z -0 .00 




90 41 + - 26 
— 0.00 






20.411H 4 


2.0 






0.673 


pn 


846 


20.68lg;f 6 




20.4ll°;g° 






20.851°-^ 


2.0 






0.842" 


pn 


1199 


21.16_° ;__ 


2.02 (10) 


20.95l°;l§ 


1.35 


(10) 


91 q 7 +0. 26 
Zl.of _ () 4 p 


9 11 +0.52 
z ' 11 -0.34 


1.06 


(9) 


0.331" 


pn 


1226 


20.8ll° ; 2 . 


0.86 (23) 


20.42l°,; 9 


0.83 


(23) 


20.81_ ; 04 o 


i - 8b _o.20 


0.78 


(22) 


0.946 


pn 


1227 


21.67l°; 2 . 




21.551°-? 






21.97l°; 32 


2.0 






0.445 


pn 


1248 


90 49+°-38 
zu -^ z -0.00 


0.97 (21) 


20 42+ ' 16 
zu -^ z -o.oo 


1.40 


(21) 


20.4ll°; 33 


1 98+ ' 23 

1 - yo -0.15 


1.01 


(20) 


0.034 


Ml M2 


1252 


20 41+° 12 

ZU.11_ go 


1.34 (16) 


20 41 +°'°° 

—0.00 


2.39 


(16) 


20 41+° 19 

ZU.11_q go 


9 9n +0.20 
z.zu_ lg 


1.24 


(16) 


0.590 


pn 


1264 


20 42+ - 28 
zu -^ z -o.oo 




90 49+ - 20 
zu.iz_ 00 






90 41 + - 81 
zu.ii_ 00 


2.0 






1.694 


pn 


1265 


99 9n+° 15 
zz.zu_ 16 




22 11+° 16 

— 0.18 






22.40l°;i? 


2.0 






0.205" 


pn 



a : spectroscopic redshift 
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